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Foreword 


These  two  mathematical  discussions  of  traffic  analyses  tend  to 
confirm  the  theory  that  if  we  are  properly  to  understand  traffic 
we  must  learn  more  about  its  behavior;  and  that  if  we  expect 
to  control  traffic  efficiently  and  provide  suitable  means  for  con- 
trol, we  must  be  able  to  uniformly  predict  traffic  patterns. 

To  know  what  to  expect  will  provide  a  guide  to  a  rule  of 
action  to  fit  the  problem. 

The  authors,  Mr.  Gerlough  and  Mr.  Schuhl,  live  in  distantly 
separated  parts  of  the  world.  It  is  safe  to  assume  they  were 
motivated  in  preparing  their  articles  by  their  daily  interest  and 
activities  in  traffic. 

They  are  highly  competent  and  have  that  very  valuable  ad- 
vantage of  experience  so  essential  in  translating  theory  to 
practice. 

Continued  efforts  in  this  direction  tend  to  provide  an  ac- 
cumulation of  basic  knowledge  of  traffic  performance  and  in 
time  to  evolve  new  concepts  of  our  traffic  problem  with  im- 
proved ideas  in  planning  and  control. 

Trusting  that  this  brochure  will  be  some  encouragement  to 
further  research,  the  Eno  Foundation  welcomes  the  opportunity 
to  publish  it  and  wishes  to  express  its  appreciation  to  the 
authors. 

Eno  Foundation 
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Use  of  Poisson  Distribution 
in  Highway  Traffic 

Daniel  L,  Gerlough 

Mr.  Gerlough,  a  Registered  Professional  Electrical 
Engineer  in  the  state  of  California,  is  an  engineer 
with  the  Institute  of  Transportation  and  Traffic 
Engineering  at  the  University  of  California  at  Los 
Angeles.  For  ten  years  prior  to  his  joining  the  In- 
stitute at  Los  Angeles  in  i()48,  he  served  as  an  engineer 
for  six  corporations  in  Texas  and  California.  Mr. 
Gerlough  is  a  member  of  the  Institute  of  Radio 
Engineers,  of  the  American  Institute  of  Electrical 
Engineers,  and  of  the  Institute  of  Traffic  Engineers. 

Applications  of  probability  and  statistics  to  engineering  prob- 
lems have  increased  markedly  during  the  past  decade.  These 
techniques  are  particularly  helpful  in  the  study  of  vehicular 
traffic.  A  monograph^  published  recently  by  the  Eno  Founda- 
tion presents  applications  of  the  broad  field  of  statistics  to 
traffic  problems.  The  present  paper  discusses  in  detail  one  par- 
ticular probability  distribution  which  has  been  found  useful  in 
the  treatment  of  several  types  of  problems  encountered  by  the 
traffic  engineer. 

The    Poisson    distribution    is    a    mathematical    relationship 
which  finds  such  applications  as: 

Testing  the  randomness  of  a  given  set  of  data 
Fitting  of  empirical  data  to  a  theoretical  curve 
Prediction  of  certain  phenomena  from  basic  data 

1  Greenshields,  Bruce  D.,  and  Weida,  Frank  M.,  Statistics  with  Applications 
to  Highway  Traffic  Analyses,  The  Eno  Foundation,  1952. 
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Some  of  the  uses  open  to  the  traffic  engineer  include: 

Analysis  of  arrival  rates  at  a  given  point 
Studies  of  vehicle  spacing  (gaps) 

Determination  of  the  probability  of  finding  a  vacant  park- 
ing space 
Studies  of  certain  accidents 
Design  of  left-turn  pockets 

The  application  of  the  Poisson  distribution  to  traffic  prob- 
lems is  not  new.  Certain  applications  were  discussed  by  Kinzer^ 
in  1933,  Adams^  in  1936,  and  Greenshields*  in  1947.  It  is  the 
purpose  here  to  provide  a  more  extensive  treatment  in  a  form 
that  will  serve  as  a  reference  guide  for  both  students  and 
practitioners  of  traffic  engineering.  In  attempting  to  fulfill  this 
objective  the  following  material  has  been  included: 

Examples  which  cover  typical  cases  and  at  the  same  time  indicate 
the  variety  of  applicability. 

Detailed  descriptions  of  the  steps  to  be  followed  in  applying  the 
Poisson  distribution  and  in  testing  it  by  the  -^  (chi-square)  test. 

Charts  for  simplified  computations  as  well  as  references  to  stand- 
ard tables  and  other  computational  aids. 

Derivation  of  the  Poisson  distribution  (in  sections  B  to  E  of  the 
Appendix)  starting  from  the  fundamental  concepts  of  permutations 
and  combinations. 

Bibliographical  references  to  more  advanced  applications  in 
which  the  Poisson  distribution  is  combined  with  other  relationships. 

2  Kinzer,  John  P.,  Application  of  the  Theory  of  Probability  to  Problems  of 
Highway  Traffic,  thesis  submitted  in  partial  satisfaction  of  requirements  for 
degree  of  B.C.E.,  Polytechnic  Institute  of  Brooklyn,  June  1,  1933. 

Abstracted  in:  Proceedings,  Institute  of  Traffic  Engineers,  vol.  5,  1934,  pp.  118- 
124. 

3  Adams,  William  F.,  "Road  TraflSc  Considered  as  a  Random  Series,"  Journal 
Institution  of  Civil  Engineers,  vol.  4,  Nov.,  1936,  pp.  121-130-f . 

4  Greenshields,  Bruce  D.;  Shapiro,  Donald;  Ericksen,  Elroy  L.,  Traffic  Per- 
formance at  Urban  Street  Intersections,  Technical  Report  No.  1,  Yale  Bureau 
of  Highway  Traffic,  1947. 
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Nature  of  Poisson  Distribution 


When  events  of  a  given  group  occur  in  discrete  degrees  (heads 
or  tails,  i  to  6  on  the  face  of  a  die,  etc.)  the  probability  of  occur- 
rence of  a  particular  event  in  a  specified  number  of  trials  may 
be  described  by  the  Bernoulli  or  binomial  distribution  (See 
Appendix,  Sections  B  and  C). 

As  an  example,  let  us  consider  an  experiment  consisting 
of  five  successive  drawings  of  a  ball  from  an  urn  contain- 
ing uniformly  mixed  black  and  white  balls,  with  the  drawn 
ball  being  returned  to  the  urn  after  each  drawing.  Of  the  five 
drawings  which  comprise  a  single  experiment,  let  x  be  the 
number  which  produced  black  balls;  thus  x  can  equal  o,  1,2,  3, 
4,  or  5.  Let  P{x)  be  the  probability  that  in  a  given  experiment 
the  number  of  black  balls  drawn  would  be  exactly  x.  If  ^  is  the 
probability  that  a  particular  drawing  will  yield  a  black  ball, 
and  q(q  =  1 — p)  is  the  probability  that  a  particular  drawing 
yields  a  white  ball,  then 

where  Ct  is  the  number  of  combinations  of  5  things  taken  x  at 
a  time.  It  can  be  seen  that  in  the  above  example  p^  the  prob- 
ability that  a  single  drawing  yields  a  black  ball,  is  equal  to 
the  percentage  of  black  balls  in  the  urn.  With  the  ball  being 
replaced  after  each  drawing  p  will  remain  constant  from  draw- 
ing to  drawing. 

The  experiment  just  described  is  an  example  of  the 
"Bernoulli"  or  "binomial"  distribution  in  which  P{x),  the  prob- 
ability of  exactly  x  successes  out  of  n  trials  of  an  event  where 
the  probability  of  success  remains  constant  from  event  to  event, 
is  given  by 

P(x)   =   Qp^q^-^ 

If  the  number  of  items  in  the  sample  n^  becomes  very  large 
while  the  product  pn  =  m  is  3.  finite  constant,  the  binomial 
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distribution  approaches  the  Poisson  distribution  as  a  limit.* 


X! 


where  e  is  the  Napierian  base  of  logarithms. 

The  derivation  of  the  Poisson  distribution  as  a  limiting  case 
of  the  binomial  distribution  is  given  in  Appendix  D.  (The 
Poisson  distribution  can  be  derived  independently  of  the  bi- 
nomial distribution  by  more  advanced  concepts.  See  Fry.^) 

The  mathematical  conditions  of  an  infinite  number  of  trials 
and  an  infinitesimal  probability  are  never  achieved  in  practical 
problems.  Nevertheless,  the  Poisson  distribution  is  useful  as 
approximating  the  binomial  under  appropriate  conditions. 

For  such  practical  purposes,  then,  the  Poisson  distribution 
may  be  stated  as  follows: 

If  in  a  given  experiment  the  number  of  opportunities  for  an 
event  to  occur  is  large  (e.g.,  n  ^  50) 

and 

If  the  number  of  times  a  particular  event  occurs  is  small 
(e.g.,  P^  0.1) 

and 

If  the  averas^e  number  of  times  the  event  occurs  has  a  finite 
value,  m,  (m  =  np), 

Then      P{x)  =  ^^- 


where  x  =  o,  1,  2,  .  .  . 
In  this  statement  of  the  Poisson  distribution  an  experiment 
may  consist  of  such  things  as: 

a.  Observing  the  number  of  micro-organisms  in  a  standard  sample 
of  blood,  X  representing  the  number  of  micro-organisms  in  any 
one  sample,  and  n  representing  the  number  of  samples  studied. 

b.  Observing  the  number  of  alpha  particles  emitted  during  each 

*  This  implies  that  the  probability  of  occurrence,  p^  becomes  very  small. 
5  Fry,  Thornton  C,  Probability  and  Its  Engineering  Uses,  Van  Nostrand,  1928, 
pp.  220-227. 
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successive  interval  of  t  seconds.  The  number  of  such  intervals 
will  be  n,  and  Xi,  X2  .  .  .  x„  will  be  the  number  of  particles  dur- 
ing the  1st,  2nd,  .  .  .  nth  intervals. 
c.    Observing  the  number  of  blowholes  in  each  of  n  castings,   x 
representing  the  number  of  holes  in  any  one  casting. 

Examples  of  Poisson  Distributions 

The  first  record  of  the  use  of  the  Poisson  distribution  to  treat 
populations  having  the  properties  listed  is  attributed  to 
Bortkewitsch  who  studied  the  frequency  of  death  due  to  the 
"kick  of  a  horse"  among  the  members  of  ten  Prussian  Cavalry 
Corps  during  a  period  of  20  years.  Some  of  the  earliest  engineer- 
ing problems  treated  by  the  Poisson  distribution  were  telephone 
problems.  The  following  example  is  based  on  such  data:^ 

Example  1 
Connections  to  Wrong  Number 

Number  of  wrong  Number  of  periods  exhibiting 

connections  per  the  given  number  of  wrong 

period  connections  (observed  frequency) 

0  o 

1  o 

2  1 

3  5 

4  11 

5  14 

6  22 

7  43 

8  31 

9  40 

10  35 

11  20 

12  18 

13  12 

14  7 

15  6 

16  2 

6  Thorndike,    Frances,    "Applications    of    Poisson's    Probability    Summation," 
Bell  System  Tedmical  Journal,  vol.  5,  no.  4,  Oct.,  1926,  pp.  604-624. 
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The  fitting  of  the  Poisson  distribution  to  the  experimental  data 
may  be  carried  out  in  tabular  form  as  follows: 

Col.  1  Col.  2  Col.  3  Col.  4  Col.  5 


X  = 

Number  of 

wrong 

Observed 

Total 

connections 

frequency 

wrong 

P{x)  = 

Theoretical 

per  period 

(periods) 

connections 

Probability 

freqiiejicy 

0 

0 

0 

0.000 

0.0 

1 

0 

0 

O.OOl 

0.3 

2 

1 

2 

0.006 

1.6 

3 

5 

15 

0.018 

4.8 

4 

11 

44 

0.039 

10.4 

5 

14 

70 

0.068 

18.2 

6 

22 

132 

0-099 

26.4 

7 

43 

301 

0.124 

.     33-1 

8 

31 

248 

0.135 

36.0 

9 

40 

360 

0.132 

35-2 

10 

35 

350 

0.115 

30-7 

11 

20 

220 

0.091 

24-3 

12 

18 

216 

0.067 

17-9 

13 

12 

156 

0.045 

12.0 

14 

7 

98 

0.028 

7-5 

15 

6 

90 

0.016 

4-3 

16 

2 

32 

0.009 

2.4 

>i6 

0 
267 

0 

0.007 
1.000 

1-9 

Total 

2334 

267.0 

The  entries  in  Col.  3  are  the  products  of  the  corresponding  entries 
in  Col.  1  and  Col.  2. 

m  =  average  number  of  wrong  connections  per  period  = 

2334 


P{x) 


^      =  S.742 
267  '^ 

(8.742)^^-8- 


In  the  last  row  of  the  table  (for  more  than  16  wrong  connec- 
tions per  period)  the  values  of  probability  and  theoretical  fre- 
quency are  necessary  to  make  the  totals  balance.  These  values 
represent  the  summation  from  1 7  to  infinity. 
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Methods  of  evaluating  e-'',  xl,  and  P{x)  are  discussed  later  in 
the  paper. 

Calculated  Frequency  =  (Total  periods  observed)  P(x) 

=  267  P{x) 

It  will  be  seen  that  there  is  a  high  degree  of  agreement  be- 
tween the  observed  and  calculated  frequencies. 

Note:  Column  5  shows  the  theoretical  or  calculated  fre- 
quencies. When  these  frequencies  are  calculated,  fractional 
(decimal)  values  often  result.  Here  these  values  have  been 
rounded-off  to  the  nearest  0.1.  The  observed  frequencies  will, 
of  course,  always  be  integral  numbers. 

Following  the  pioneer  work  in  the  field  of  telephone  appli- 
cations, the  Poisson  distribution  was  gradually  applied  to  other 
engineering  problems.  The  following  example  adapted  from 
Grant^  shows  an  application  to  the  occurrence  of  excessive 
rainfall : 

Example  2 


Rainstorms 

Number 

"j 

of  storms 

per 

0 bserved 

Theoretical 

station 

number  of 

Total 

number  of 

per  year 

occurrences 

storms 

Probability 

occurrences 

0 

102 

0 

0.301 

99-3 

1 

114 

114 

0.361 

119.1 

2 

74 

148 

0.217 

71.6 

3 

28 

84 

0.087 

28.7 

4 

10 

40 

0.026 

8.6 

5 

2 

10 

0.006 

2.0 

>5 

0 

330 

0 

396 

0.002 

0.7 

Total 

1.000 

330.0 

m  _  396  _ 

1.20 

e-"^  = 

0.301 

330 

'  e-""  =  99-' 

330 

The  first  published  examples  of  the  Poisson  distribution  as 

7  Grant,    Eugene    L.,    "Rainfall    Intensities    and    Frequencies,"    Transactions, 
American  Society  of  Civil  Engineers,  vol.  103,  1938,  pp.  384-388. 


8 


POISSON    DISTRIBUTION 


applied  to  traffic  data  are  presented  by  Adams.^  The  following 
is  one  of  Adams'  examples: 

Example  3 

Rate  of  Arrival  (Vera  St.) 
(Number  of  vehicles  arriving  per  10  second  interval) 


No. 

vehicles 

per  10 

' 

sec. 

period 

Observed 

Total 

Probability 

Theoretical 

X 

frequency 

vehicles 

P(x) 

frequency 

0 

1 

94 
63 

0 
63 

0-539 
0-333 

97-0 
59-9 

2 

21 

42 

0.103 

18.5 

3 

2 

6 

0.021 

3.8 

>3 

0 

180 

0 

111 

0.004 
1.000 

0.8 

Total 

180.0 

m 

total  vehicles        1 1 1 
total  periods         180 

=  0.617;  ^~-^"  = 

•539 

Pi.\  _  (-^ 

'7)'   ,-.«7   , 

(•539)  (-617)" 

Calculated  frequency  =180  P(x) 

Note:  Since  there  were  1 1 1  vehicles  in  180  ten-second  periods, 
the  hourly  volume  was  222. 


Testing  Goodness  of  Fit  (^^  Test) 

In  each  of  the  foregoing  examples  it  has  been  postulated  that  a 
Poisson  distribution  having  a  parameter  m  whose  value  has 
been  computed  from  the  observed  data  describes  the  popula- 
tion that  has  been  sampled.  The  observed  distribution  con- 
stitutes this  sample.  By  inspection  there  is  apparent  agreement 
between  the  observed  distribution  (sample)  and  the  theoretical 
distribution.  The  inference  is  then  made  that  the  postulated 
theoretical  (Poisson)  distribution  is  in  fact  the  true  population 

8  Adams,  op.  cit. 
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distribution.  This  inference  is  based,  however,  solely  on  inspec- 
tion; a  more  rigorous  basis  for  reaching  such  a  conclusion  is 
desired.  One  of  several  statistical  tests  of  significance  may  be 
used  for  this  purpose;  the  chi  square  (x^)  test  is  appropriate  to 
the  present  application.  This  test,  which  is  described  in  Ap- 
pendix G  and  illustrated  below,  provides  for  one  of  two 
decisions: 

1.  It  is  not  very  likely  that  the  true  distribution  (of  which  the  ob- 
served data  constitute  a  sample)  is  in  fact  identical  with  the 
postulated  distribution. 

2.  The  true  distribution  (of  which  the  observed  data  constitute  a 
sample)  could  be  identical  with  the  postulated  distribution. 

It  can  be  seen  that  either  decision  can  be  erroneously  made. 
Decision  i  can  be  wrong  if  in  fact  the  postulated  distribution 
is  the  true  distribution.  On  the  other  hand.  Decision  2  can  be 
wrong  if  the  true  distribution  is  in  fact  different  from  the  pos- 
tulated distribution.  Statistical  tests  of  significance  allow  for 
specifying  the  probability  of  making  either  of  these  types  of 
error.  Usually  the  probability  of  making  the  first  type  of  error 
(incorrectly  rejecting  the  postulated  distribution  when  in  fact 
it  is  identical  with  the  true  distribution)  is  specified  and  no 
statement  is  made  with  regard  to  the  second  type  of  error.  The 
specification  of  the  first  type  of  error  is  expressed  as  a  "signifi- 
cance level."  Common  significance  levels  are  0.01,  0.05,  and 
0.10.  Thus,  when  a  test  is  made  at  the  0.05  (5%)  level,  the 
engineer  takes  the  chance  that  5%  of  the  rejected  postulated 
distributions  are  in  fact  identical  with  the  corresponding  true 
distributions.  For  a  complete  discussion  of  the  theory  under- 
lying this  and  other  statistical  tests  of  significance,  the  reader 
is  referred  to  any  standard  text  on  statistics. 

The  technique  of  performing  the  x^  significance  test  is  illus- 
trated in  Example  4.  In  the  table  of  Example  4  the  observed 
frequency  of  wrong  connections  is  shown  along  with  the  pos- 
tulated theoretical  frequency  distribution  (computed  in  Ex- 
ample 1).  The  problem  is  to  decide  whether  the  observed  data 
can   be  construed   as   a  sample   coming  from   the   postulated 
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theoretical  distribution.  The  computing  techniques  described 
in  Appendix  G  yield  a  chi  square  value  of  7.6  with  1 1  degrees  of 
freedom.  (It  should  be  noted  in  the  table  of  Example  4  that  the 
frequencies  of  3  or  fewer  wrong  connections  have  been  com- 
bined as  have  the  frequencies  of  15  or  greater.  This  is  done  to 
meet  the  requirement  of  the  x^  test  that  the  theoretical  fre- 
quency be  at  least  5  in  any  group.  Thus,  this  combination  re- 
sults in  13  groups,  with  13  —  2^11  degrees  of  freedom.) 

Example  4 
X^  Test  of  Connections  to  Wrong  Number  (Data  from  Example  1) 
Number  of 


wrong  con- 

Postulated 

nections  per 

Observed 

theoretical 

period 

frequency 

frequent 

cy 

P 

X 

f 

F 

F 

^3 

6 

6.7 

5-4 

4 

11 

10.4 

11.6 

5 

14 

18.2 

10.8 

6 

22 

26.4 

18.3 

7 

43 

33-1 

55-9 

8 

31 

36.0 

26.7 

9 

40 

35-2 

45-5 

10 

35 

30.7 

39-9 

11 

20 

24-3 

16.5 

12 

18 

17-9 

18.1 

13 

12 

12.0 

12.0 

14 

7 

7-5 

6.5 

—  15 

8 

267 

8.6 

7-4 

fOTAL 

267.0 

274.6 

X^  = 

274.6  - 

-  267.0  = 

7.6 

From  statistical  tables  or  from  Figure  3  the  value  of  x^  at 
the  0.05  (5%)  significance  level  for  v  =  1 1  is  found  to  be  19.7.* 

*  In  using  tables  of  p^2  care  should  be  exercised  to  note  the  manner  in  which 
the  table  is  entered  with  the  significance  level.  If  the  table  is  so  constructed 
that  for  a  given  number  of  degrees  of  freedom  the  value  of  ^2  increases  with 
decreasing  percentiles  (probabilities),  the  table  is  entered  with  the  percentile 
corresponding  to  the  significance  level.  If  the  table  is  such  that  the  value  of  -^^ 
increases  with  increasing  percentiles,  the  table  is  entered  with  the  significance 
level  subtracted  from  i. 


POISSON    AND    TRAFFIC  1  1 

This  is  known  as  the  critical  value  for  the  significance  test  at 
the  5%  level.  Since  the  computed  value  of  x^  (7-6)  is  less  than 
the  critical  value,  decision  2  is  accepted.  (The  acceptance  of 
decision  2  is  rigorously  stated,  "There  is  no  evidence  to  indicate 
that  the  true  theoretical  distribution  differs  from  the  postulated 
distribution.")  Had  the  computed  value  of  x^  exceeded  the 
critical  value,  decision  1  would  have  been  accepted. 

Once  the  decision  has  been  made  that  a  particular  postulated 
theoretical  distribution  appropriately  represents  the  population 
from  which  the  observed  data  came,  the  theoretical  distribution 
can  be  used  in  place  of  the  observed  distribution  for  engineer- 
ing analysis  and  action,  since  it  will  be  free  from  the  random 
variations  present  in  the  observed  distribution. 

Interpretation  of  Examples  1  to  3 

When  the  engineer  gathers  field  data  he  is  usually  interested  in 
predicting  performance  for  design  or  administrative  purposes. 
Fitting  a  theoretical  curve  to  empirical  data  is  only  one  step  in 
the  overall  objective.  The  Poisson  distribution  is  of  value  only 
if  it  permits  useful  conclusions  to  be  drawn.  Below  are  some 
conclusions  which  can  be  drawn  from  each  of  examples  1  to  3. 

Example  i: 

a.  Wrong  connections  follow  a  Poisson  distribution  with  the 
parameter  m  =  8.742. 

b.  Probabilities  computed  from  the  Poisson  distribution  may  be 
used  to  predict  the  occurrence  of  wrong  connections  and  thus 
may  serve  as  a  basis  for  future  design  considerations.  Future 
design  might,  for  instance,  provide  additional  facilities  to  com- 
pensate for  the  capacity  consumed  by  the  wrong  connections. 
A  possible  alternative  might  be  the  provision  of  some  method 
of  preventing  wrong  connections. 

c.  The  effectiveness  of  engineering  changes  may  be  evaluated  in 
terms  of  a  significant  change  in  Poisson  parameter. 

Example  2: 

2L.  The  occurrence  of  rainstorms  in  the  population  considered 
follows  the  Poisson  distribution  with  the  parameter  m  =  1.20. 
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b.  The  probabilities  may  serve  as  a  basis  for  future  design. 
Consider  for  instance  the  following: 

Let  Di      =  the  damage  caused  if  one  storm  occurs  during  a 
year, 
D2      =  the  damage  caused  if  two  storms  occur  during  a 
year, 

Dn      =  the  damage  caused  if  n  storms  occur  during  a 

year, 
P(i)   =  the  probability  of  one  storm  during  a  year, 
P(2)   =  the  probability  of  two  storms  during  a  year, 

P(n)  =  the  probability  of  n  storms  during  a  year. 
Then  the  expected  or  average  damage  per  year  will  be 
Da,.  =  £>i  P(i)  +  D2  P{2)  + +  D,  P{n). 

Example  3.* 

a.  At  the  location  considered  and  during  the  time  of  day  and  day 
of  the  week  considered  the  arrival  of  vehicles  follows  the 
Poisson  distribution  with  the  parameter  m  =  0.617  ^^^  i^- 
second  intervals. 

b.  The  probabilities  may  be  used  for  design  of  traffic  control 
measures.  (Examples  appearing  later  in  this  paper  indicate 
some  specific  methods  of  applying  probabilities  to  traffic  de- 
sign problems.) 

Traffic  Applications  of  the  Poisson  Distribution 

Prediction  of  Arrivals  from  Hourly  Volume 

For  conditions  in  which  the  Poisson  distribution  applies— i.e., 
under  conditions  of  "free  flow"— it  is  possible  to  compute  the 

probability  of  o,  1,  2, ,  k  vehicles  arriving  per  time 

interval  of  t  seconds  provided  the  hourly  volume,  Vj  is  known: 

t     =  length  of  time  interval  in  seconds 
V  =  hourly  volume 
n    =  number  of  intervals  (per  hour) 
3600 
^      t 
m  =  average  number  of  vehicles  per  interval 

V_  _     Vt 

3600         3600 
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Then  the  probability,  P{x),  that  x  vehicles  will  arrive  during 
any  interval  is: 

The  hourly  frequency.  Fa,,  of  intervals  containing  x  vehicles  is: 


vt 

3600 


If  the  period  under  consideration  is  different  from  an  hour,  the 
3600  in  the  first  parentheses  would  be  replaced  by  the  appro- 
priate length  of  time  in  seconds.  The  value  V,  however,  would 
still  be  the  hourly  volume.  Examples  5  and  6  illustrate  this. 

Example  5 
Prediction  of  Arrivals  (Low  Volume) 
V   =  Hourly  volume  =  37 
t     =  Length  of  interval  =  30  seconds 

_  _    Vt    _  37  •  30 


3600  3600 

0-735 

3600  _  3600 


0.308 


120 


t  30 

Fo  =  n  l^\  m°  e-^  =  (120)  (i)  (i)  (0.735)  =  88.2 

Fi  =  Fo  f  —  j  =  (88.2)  (0.308)  =  27.2 


F,  =  F.  {^\  =  &MJ^ 


308) 

- — -  =  4.2 


F  >2=  n  -  (F0  +  F1  +  F2)  =  120  -  (88.2  +  27.2  +  4.2)  =  0.4 

Note:  The  method  of  computing  F-^  from  Fq  and  F2  from  Fj  is  dis- 
cussed in  section  on  Methods  of  Computation, 
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The  following  table  compares  the  predicted  frequencies  with 
frequencies  observed  in  a  field  study  (Elliot  Avenue  West- 
bound).^ The  predicted  frequency  in  this  example  corresponds 
to  the  theoretical  frequency  of  the  previous  examples. 


Number  of 

vehicles 

arriving  per 

Observed^ 

Predicted 

interval 

frequency 

frequency 

X 

fa> 

F. 

o 

88 

88.2 

1 

27 

27.2 

2 

5 

4.2 

^3 

o 

0.4 

120 


120.0 


Example  6 

Prediction  of  Arrivals  (Moderate  Volume) 
(30-Second  Intervals  for  25  Minutes) 

V   =  418 

t     =30  seconds 

418  •  30  _ 

2^^   •   60  1500  /  •        ^  r  1  •  N 

n    =  -^ =  -^ —  =  50  (2^  minutes,  00  seconds  per  minute) 

30  30  ^     '  ^ 

e-^=  0.0308 

Fo  =  n  e~^  =  50  •  0.308  =   1.54 

etc. 


In  the  following  table  the  various  predicted  frequencies  are 
compared  with  frequencies  observed  in  the  field  (Redondo 
Beach  Blvd.  Eastbound). 

9  Based  on  field  data  supplied  through  the  courtesy  of  the  Traffic  and  Lighting 
Division,  Los  Angeles  County  Road  Department. 
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Number  of 

vehicles 

arriving  per 

Observed 

Predicted 

interval 

X 

frequency 

u 

frequency 
F. 

0 

0 

1-5 

1 

9 

5-4 

2 

6 

9-3 

3 

9 

10.8 

4 
5 
6 

11 

9 
5 

9-4 
6.6 

3.8 

7 
^8 

1 
0 

50 

1-9 

1-3 

50.0 

The  following 

table  shows  the 

X"  analysis  of  the  data  abo^ 

Number 

vehicles 

Predicted 

arriving  per 
interval 

X 

Observed 
frequency 

u 

(theoretical) 
frequency 

F. 

Fa, 

^1 

9 

6.9 

11.7 

2 

6 

9-3 

3-9 

3 

9 

10.8 

7-5 

4 

5 

^6 

11 

9 
6 

9-4 
6.6 
7.0 

12.9 

12.3 

5-1 

50 

50.0 

53-4 

M- 

=  53-4  - 

■  50.0  = 

3-4 

V 

=  6-2 

=  4 

Xl)5 


The  agreement 
ceptable  at  the 


=  9.49  (The  symbol  x^os  is  used  to  denote  the  critical 
value  of  x^  at  the  0.05  (5%)  significance  level, 
between  the  predicted  and  observed  data  is  ac- 
5%  significance  level. 


Testing  Goodness  of  Fit  with  a  Poisson  Population 

In  some  locations   traffic  parameters   will   vary  with   time 


i6 
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throughout  the  day  (or  week  or  year)  while  at  other  locations 
the  variability  may  appear  to  be  entirely  random.  One  way 
in  which  the  traffic  engineer  may  test  the  randomness  of  his 
data  is  to  assume  that  it  follows  the  Poisson  distribution,  and 
then  to  test  this  assumption  by  the  x^  test.  If  the  x^  test  fails, 
there  is  strong  evidence  that  the  data  are  non-random.  If  the 
results  of  the  x^  test  are  acceptable,  the  data  may  be  assumed 
to  be  random  for  the  purposes  of  the  traffic  engineer.  It  should 
be  noted,  however,  that  an  acceptable  x^  test  is  less  conclusive 
than  one  which  fails.  The  rigorous  interpretation  of  an  ac- 
ceptable test  should  be,  "There  is  no  evidence  to  indicate  non- 
randomness  of  the  data." 

Example  7  is  an  analysis  of  right  turns  during  300  3-minute 
intervals  distributed  throughout  various  hours  of  the  day  and 
various  days  of  the  week.  While  casual  inspection  of  the  ob- 
served data  might  lead  one  to  suspect  randomness,  the  test 
shows  strong  evidence  of  non-randomness. 


Example  7 

Right  Turns 

(Westwood  at  Pico) 

Number  of  cars 

making  right 

Obseii~ved 

Theoretical 

turns  per  ymin. 

frequency 

frequency 

P 

interval 

f 

F 

F 

0 

14 

6.1 

32.1 

1 

30 

23-7 

38.0 

2 

36 

46.2 

28.1 

3 

68 

59-9 

77.2 

4 

43 

58.3 

31-7 

5 

43 

45-4 

40.7 

6 

30 

294 

30.6 

1 

14 

16.4 

12.0 

8 

10 

8.0 

12.5 

9 

6" 

3-4 

10 

4 

1-3 

11 

1 

-12 

0.6 

.6.6 

21.8 

12 

1 

0.3 

^13 

0 

1.0 

300 


300.0 


524-7 


(The  brackets  indicate  grouping  for  the  x^  test.) 


I 


POISSON    AND   TRAFFIC 

ii68 

300         ^    ^^ 

F 

-  n  =  324.7  -  300.0  =  24.7 

1;  =  10  —  2  =  8 

X^05    =    15-5 

17 

0.0203 


Since  24.7  >  15.5,  the  fit  is  not  acceptable  at  the  5%  signifi- 
cance level. 

Some  traffic  phenomena  may  be  random  when  observed  for  an 
interval  of  one  length  but  non-random  when  observed  with  an 
interval  of  a  different  length.  Example  8a  shows  arrival  rates 
for  lo-second  intervals.  Here  the  failure  of  the  x^  test  (non- 
acceptability  of  fit)  indicates  non-randomness.  Example  8b  con- 
siders the  same  arrivals  with  30-second  observation  intervals; 
the  assumption  of  randomness  acceptable  at  the  5%  level. 


Example 

8a 

VAL  RaTE- 

-io-Second  Intervals 

(Durfee  Avenue, 

Northb 

Jumber  of 

Observed 

Theoretical 

cars  per 

frequency'^^ 

frequency 

f 

interval 

f 

F 

F 

0 

139 

129.6 

149.1 

1 

128 

132.4 

123.7 

2 

55 

67.7 

44-7 

3 

25 

23.1 

27.1 

4 

10 

5-9'" 

5 

3^3 

1.2  I7. 2 

23.5 

^6 

oj 

0.1 

360  360.0  368.1 

(As  before,  the  brackets  indicate  grouping  for  the  ^^  test) 

368 
m  =  ^V~  =   1.022  e~^  =  o.ribo 

360 

10  See  footnote  9. 
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^~  —  n  =  368.1   —  360.0  =  8.1 

i;   =    5    -    0    =   3 

x'o5  =  7-8i 


The  fit  is,  therefore,  not  acceptable  at  the  5%  level. 

Example  8b 
Arrival  Rate— 30-SECOND  Intervals  (Durfee  Avenue,  Northbound) 


Number  of 

Obsented 

1  heoretical 

cars  per 

jrequency'^'^ 

frequency 

f 

interval 

f 

F 

F 

0 

9 

5.6 

14-5 

1 

16 

17.2 

14-9 

2 

30 

26.3 

34-2 

3 

22 

26.9 

18.0 

4 

19 

20.6 

17-5 

5 

10 

12.6 

7-9 

6 

5 

6-51 

1 
8 

1 
3 

-14 

2.8 
1.1 

10.8 

i8.i 

^9 

1 

0.4 

120 

120.0 

125.1 

m 

_  368  _ 

120 

=  3-f 

(67 

g-m    — 

0.047 

I 


r- 


—  n  =   125. 1   —   120.0  =  5.1 


X^05 


7-2=5 
I  I.I 


This  fit  is  acceptable  at  the  5%  significance  level. 

Probability  of  Fiiiding  a  Vacant  Parking  Space 

Certain  analyses  of  parking  may  be  treated  by  the  Poisson 
distribution.  Example  9  shows  the  results  of  a  study  of  4  block 


10  See  footnote  9. 
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faces  containing  48  one-hour  parking  spaces.  (Loading  zones 
and  short-time  parking  spaces  were  not  included).  Observations 
were  taken  during  the  hours  of  2  P.M.  to  4  P.M.  on  the  days 
Monday  through  Friday.  The  observation  period  was  divided 
into  5-minute  intervals.  During  each  5-minute  interval  one 
observation  was  made;  the  exact  time  within  each  interval  was 
randomized.  Each  observation  consisted  of  an  instantaneous 
count  of  the  number  of  vacant  spaces.  Since  the  x^  test  shows 
good  agreement  between  the  observed  and  theoretical  data, 
there  is  no  ground  to  doubt  that  the  distribution  of  vacant 
parking  spaces  follows  the  Poisson  distribution.  Assuming  the 
Poisson  distribution  to  hold,  it  is  possible  to  compute  the 
probability  of  finding  a  vacant  one-hour  parking  space  within 
these  four  block  faces. 

Example  9 

Vacant  Parking  Spaces  (Westwood  Blvd.) 
Number  of 


t 
F 

33-5 

44-9 

14.7 

15.6 
15.8 


vacant  i 

-hour 

parking 

Observed 

Theoretical 

spaces  per 
observation 

frequency 
f 

frequency 
F 

0 

29 

25.1 

1 

42 

39-3 

2 

21 

30.1 

3 

16 

16.4 

4 

7 

6.4 

5 

6 

^7 

2 

3 
0 

'12 

2.0 
0.5 
0.2 

■9-1 

120 


120.0 


124.5 


I 


F 


m  = 


V  = 
2     _ 


x:o5 


120 


1.567 


n  =  124.5  "  120.0  =  4.5 
5-2  =  3 


7.81 


0.209 


20 


POISSON    DISTRIBUTION 


The  probability  of  finding  one   or  more  vacant  parking  spaces, 
P  (^  i),  is 

P  (^  i)  =  1  —  P{0)  =  1  —  e-'"  =  1  —  0.209  =  0.791 

Analysis  and  Comparison  of  Accident  Data 

Several  writers  have  discussed  the  use  of  the  Poisson  dis- 
tribution in  accident  analysis.  The  example  which  follows  is 
taken  from  an  unpublished  paper  by  Belmont.^^ 

Ninety-three  road  sections  each  having  average  traffic  less 
than  8,000  vehicles  per  day  and  average  usage  of  less  than 
40,000  vehicle  miles  per  day  were  combined  to  give  45  com- 
posite roads  each  carrying  35,000  to  40,000  vehicle-miles  per 
day.  The  accident  records  of  these  roads  were  then  compared 
usingr  the  Poisson  distribution  as  follows: 

Example  10 
Single-car  Accidents 


Number  of 

single-car 

Number  of 

Theoretical 

accidents 

roads 

number  of 

ip^O 
X 

observed 

u 

roads 

Fa, 

F, 

0 

i8 

14-5 

22.3 

1 

14 

16.4 

12.0 

2 

1 

9-3 

3 

4 

3-5 

4 

1 

1.0 

5 

0 

13 

0.2 

-14.1              12.0 

6 

0 

0.04 

7 
^8 

1 
0 

0.007 
0.001 

45 

45 

46.3 

m  = 

51 
45 

=  I.I 

33 

^-m    _    o_222 

45^-^  =  ' 

14 


•5 


11  Belmont,  D.  M.,  The  Use  of  the  Poisson  Distribution  in  the  Study  of  Motor 
Vehicle  Accidents,  Technical  Memorandum  No.  B-7,  Institute  of  Transportation 
and  TrafiSc  Engineering,  University  of  California,  June  20,  1952. 
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I 


46.3  -  45.0  =  1.3 


i;   =    3    -   2 


X!)5   =    3-84 


Although  the  fit  is  acceptable  at  the  5%  level,  one  might  be 
suspicious  of  the  road  section  having  7  accidents  when  the 


probability  of  this  occurrence  is 


0.007 

45 


or  approximately  2  X 


io~^.  Eliminating  this  one  road  results  in  the  following: 

Example  10b 
Single-car  Accidents 


Number  of 

single-car 

Number  of 

Theoretical 

accidents 

roads 

number  of 

1950 

observed 

roads 

f 

X 

fee 

Fa, 

J^. 

0 

18 

16.2 

20.0 

1 

14 

16.2 

12.1 

2 

7" 

8.r 

3 

4 
^5 

^       4 

0 

44 

2.7 
0.7 
0.1 

.11.6 

12.4 

44.0 

44-5 

m  = 

44 
^^  =  1. 000 

g-m    ^    Q  2 

68 

44^"^  =  I 

I 


f^ 
F. 


44 


44.5  -  44.0  =  0.5 


y  =  3  _  2  =  I 


X^05    =    3-84 

Since  there  is  no  evidence  (at  the  5%  level  of  significance)  that 
there  is  a  difference  between  the  observed  and  the  theoretical 
data,  one  may  reasonably  assume  that  each  of  these  44  roads 
has  the  same  accident  potential.  Thus,  among  road  sections 
from  this  group,  any  section  is  as  likely  as  any  other  to  have 
from  zero  to  four  accidents  per  year;  a  section  which  has  four 
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accidents  is  no  more  dangerous  than  one  havins^  no  accidents, 
the  occurrences  being  simply  the  operation  of  chance. 

When  only  two  roads  or  intersections  are  to  be  compared  as 
to  accident  potential,  the  following  technique^-  may  be  used: 

a.  Compute  u^ 

,  Xi    —    X2    —    I 

where  u 


X\  =  larger  of  observed  values 

xo  =  smaller  of  observed  values 

Compare  the  resulting  value  of  u  with  the  appropriate  critical 
value  selected  from  table.  If  the  computed  value  of  u  is  less  than 
the  critical  value,  the  two  roads  or  intersections  may  be  con- 
sidered to  have  the  same  accident  potential. 

Critical 
Significance  value  of 

level  u 

o.oi  2.58 

0.05  1.96 

o.io  1.64 

Example  11 
Comparison  of  AccroENXs  at  Two  Intersections 

Intersection  1:  6  accidents  during  year  considered 
2:  8  accidents  during  year  considered 

a:]  =  8  A-2  =  6 

8-6-1 


V8  +  6 


=  0.267 


Since  the  computed  value  of  u  in  Example  11  is  less  than  the 
5%  significance  level,  there  is  no  evidence  of  a  difference  in 
accident  potential  between  the  two  intersections. 

12  Hald,  A.,  Statistical  Theory  with  Engineering  Applications,  John  Wiley  and 
Sons,  Inc.,  1952,  pp.  725-726. 
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Cumulative  Poisson  Distribution 

In  the  previous  discussion,  the  probability  of  exactly  o,  1,  2,  .  .  . 
items  per  trial  (or  time  interval)  has  been  computed.  In  many 
problems  it  is  desirable  to  compute  such  values  as  the  prob- 
ability that  the  number  of  items  (vehicles  arriving)  per  trial  is: 

a.  k  or  less 

b.  greater  than  k 

c.  less  than  k 

d.  k  or  greater 

These  probabilities  involve  the  cumulative  Poisson  distribution 
and  may  be  expressed  as  follows: 

P{x  <  k)  =  Probability  that  x  <  k 

x  =  o  x  =  o  '^' 

P(x  >  k)  =  I   -  P(x  <  k)  =  I   -     S     ^  e-"' 

x  =  o  ^■ 

^~  ^  m^ 
P(x  <  k)  =  P(x  <  k-i)   =      S     ^^-- 

x^o  ^■ 

'^~  ^  m* 
P{x  >  k)  =   I   -  P(x  <  k)  =   I   -     S     -^e~^ 

x  =  o  ^' 

The  following  examples  illustrate  the  use  of  the  cumulative 
Poisson  distribution. 

Example  12 

Left-turn  Control 

This  example  is  patterned  after  one  given  by  Adams.^^ 
A  single  intersection  is  to  be  controlled  by  a  fixed-time  signal 
having  a  cycle  of  55  seconds.  From  one  of  the  legs  there  is  a  left- 
turning  movement  amounting  to  175  vehicles  per  hour.  The 

13  Adams,  op.  cit. 


24  POISSON    DISTRIBUTION 

layout  o£  the  intersection  is  such  that  two  left-turning  vehicles 
per  cycle  can  be  satisfactorily  handled  without  difficulty, 
whereas  three  or  more  left-turning  vehicles  per  cycle  cause 
delays  to  other  traffic.  In  what  percent  of  the  cycles  would  such 
delays  occur? 

m  =  avar.  left  turns  per  cycle  =         ^ — —  =  2.67 
^  ^        ^  3600  ' 

P{X>3)    =    1-        S        i^Ml  ^-2.67 

x  =  o       ^■ 
=  I  —  [0.069  +  0.185  +  0.247] 

=  I  -  0.501  =  0.499 

Answer:  49.9% 

If  a  special  left-turn  phase  is  provided,  in  what  percent  of  the 
cycles  will  this  special  feature  be  unnecessary  by  virtue  of  the 
fact  that  there  is  no  left-turning  vehicle? 

P{0)   =   e-2-67  =   0.069 

Answer:  6.9% 

Example  13 
Left-turn  Storage  Pockets 

The  California  Division  of  Highways  has  applied  the  cumu- 
lative Poisson  distribution  to  the  design  of  left-turn  storage 
pockets.  The  pockets  are  designed  in  such  a  manner  that  the 
number  of  cars  desiring  to  make  a  left-turn  during  any  signal 
cycle  will  exceed  the  pocket  capacity  only  4%  of  the  time— i.e., 
the  probability  of  greater  than  a  given  number  of  cars  is  0.04. 
This  may  be  expressed: 

0-04=  I  -      S    -^    -^      e    3600 


Where:  L  =  number  of  left-turning  cars  per  hour 
t    =  length  of  signal  cycle  in  seconds 
k  =  pocket  capacity  (number  of  cars) 
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This  relationship  may  be  solved  for  k  by  accumulating  terms  in 
the  right-hand  member  until  the  equation  is  satisfied.*  A  more 
direct  method  of  computation,  however,  is  to  use  tables  or 
charts  showing  the  various  values  for  the  cumulative  Poisson 
distribution.  Such  a  chart  is  included  in  the  section  on  Methods 
of  Computation. 

The  values  obtained  are  as  follows: 

Required  Storage  (Number  of  Vehicles) 


Peak  hour 

left- 

60  second 

120  second 

turn  movements 

cycle 

cycle 

lOO 

4 

7 

200 

7 

12 

300 

9 

16 

400 

12 

20 

500 

14 

24 

600 

16 

28 

700 

18 

33 

800 

20 

36 

The  Exponential  Distribution 

The  previous  discussions  of  arrival  rate,  etc.,  have  dealt  with 
the  distribution  of  discrete  events  (arrivals  of  cars)  within  a 
given  time  interval.  The  distribution  of  gaps  (time  spacing) 
between  vehicles  is  a  continuous  variable  and  is  exponential  in 
nature.  This  distribution  may  be  derived  from  the  Poisson 
distribution  as  follows: 
As  before,  let: 

m  =  average  number  of  cars  per  time  interval 

3600 

*  The  equation  will  in  general  never  be  precisely  satisfied  because  of  discrete 
nature  of  the  Poisson  distribution.  There  will  be  two  values  of  k  which  nearly 
satisfy  the  relationship— one  too  large  and  one  too  small.  It  is  then  a  matter  of 
engineering  judgment  which  to  use. 
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Where  V  =  hourly  volume 

t  =  length  of  interval  in  seconds 
P{o)  =  Probability  of  zero  vehicles  during  an  interval  t 

„,  .         m°e~^  -^finPi  (A) 

p(o)  =  -jrj-  =  e-^  =  e     ^^^^  ^    ' 

If  there  are  no  vehicles  in  a  particular  interval  of  length  t, 

then  there  will  be  a  gap  of  at  least  t  seconds  between  the  last 

previous  vehicle  and  the  next  vehicle.  In  other  words,  P{p)  is 

also  the  probability  of  a  gap  equal  to  or  greater  than  t  seconds. 

This  may  be  expressed: 

vt 
P(g>t)  =  e    3600  (B) 

From  this  relationship  it  may  be  seen  that  (under  conditions 
of  random  flow)  the  number  of  gaps  greater  than  any  given 
value  will  be  distributed  according  to  an  exponential  curve. 
This  is  known  as  an  exponential  distribution. 

Example  14 

Distribution  of  Gaps  (Arroyo  Seco  Freeway) 
Total  gaps  =  214 
Total  time  =  1753  sec. 

214 
™  —  t  —  0.122/ 


1753 

P{g>i)  =  ^~°-'''' 

G  =  (Expected  number  of  gaps  >  /)  =  2i^e~°-'^^^^ 

Figure  1  shows  a  comparison  between  the  observed  and 

theoretical  results.^^ 

Vt 
14  When  using  equation   [B]   it  should  be  noted   that  — - — ,   the  number  of 

vehicles  arriving  during  an  interval,  and  G,  the  number  of  gaps  equal  to  or 
greater  than  some  specified  value  t,  are  discrete  variables.  That  is,  the  number 

of  vehicles  or  the  number  of  gaps  can  take  on  only  integral  values  o,  1,2, 

etc.  Probabihty,  P(x)  and  time,  t,  on  the  other  hand,  are  continuous  variables 
and  can  take  on  fractional  values.  This  is  illustrated  in  Figure  1,  where  the 
theoretical  frequencies  are  represented  by  bars  and  the  probability  curve  is 
represented  by  a  continuous  dashed  curve. 


/ 


/, 


>h 
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Example  15 
Safe  Gaps  at  School  Crossings 

(This  example  is  based  on  the  report  of  a  joint  committee  of 
the  Institute  of  Traffic  Engineers  and  the  International  Asso- 
ciation of  Chiefs  of  Police.^^) 

In  studying  the  natural  gaps  in  traffic  at  school  crossings  the 
following  assumptions  may  be  made: 

1.  The  walking  speed  of  a  child  is  3.5  ft./sec. 

2.  There  must  be  at  least  one  opportunity  of  crossing  per  minute. 
(This  implies  a  minimum  of  60  opportunities  per  hour.) 

On  the  basis  of  these  assumptions  it  is  desired  to  determine  the 
critical  volume  for  a  street  of  a  given  width.  The  critical 
volume  here  refers  to  that  volume  above  which  special  measures 
will  be  required  for  the  safety  of  the  child. 

There  are  two  approaches  to  this  problem,  and  care  must  be 
exercised  in  the  selection  of  the  appropriate  method.  One  ap- 
proach is  that  of  the  number  of  gaps  greater  than  the  time  t 
required  for  a  child  to  cross  the  street.  The  expected  number 
of  gaps  per  hour  will  be  V,  and  the  probability  (fraction)  of 
gaps  equal  to  or  greater  than  t  is  as  derived  in  equation  [B]. 
Thus,  the  expected  number  of  gaps  per  hour  which  are  equal  to 
or  greater  than  t  will  be 

vt 

Yg      3600 

The  other  approach  is  the  number  of  ^second  intervals  per 
hour  which  are  free  of  cars.  The  number  of  ^second  intervals 
per  hour  is  3600/L  The  probability  (fraction)  of  such  intervals 
free  of  cars  is  as  given  in  equation  [A].  Thus,  the  expected 
number  of  ^second  intervals  per  hour  which  are  free  of  cars 
will  be 

3^,-3600 


15  "Report    on   Warrants    for   Traffic    Officers   at   School    Intersections,"   Pro- 
ceedings, Institute  of  Traffic  Engineers,  vol.  18,  1947,  pp.  118-130. 
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When  one  considers  the  fact  that  a  long  gap  may  contain 
several  ^second  intervals  during  which  there  is  opportunity  to 
cross,  this  second  approach  seems  more  appropriate. 

Introducing  the  second  assumption  (at  least  60  opportunities  to 
cross  per  hour), 

36^,-^  =  60 
t 

where  Vc  denotes  the  critical  volume. 
From  the  first  assumption  (walking  speed  =  3.5  ft/sec), 

_    D_ 

where  D  is  the  width  of  the  street.  . 
Substituting: 

3^^^1.3:5.-365^-60 
D  ~ 

12,600     210 

Takinsr  the  natural  logarithm  of  both  sides: 


-VrD        ,     D 

=  In 


12,600  210 


D  210 


Making  use  of  the  relationship 

losf  X  loff  X 

In  X  =  r^  =  — - — 
log  e         0.4343 
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gives 

12,600       I       ,        D 

Vc   = jz log  

D    0.4343    ^210 

or 

-.^  29,000    /,  r>  1  X 

Vc  = ^^  (log  D  -  log  210) 

or 

29,000  .  ,      „. 

Vc  =     ^^      (2.322  -  log  D) 

This  equation  has  been  adopted  by  the  joint  committee  of  the 
Institute  of  Traffic  Engineers  and  the  International  Association 
of  Chiefs  of  Police.  Solution  of  the  equation  gives  the  following 
values 

Width  of  street  (ft.):  25  50  75 

Critical  volume  (v.p.h.):    1072  361  173 

Miscellaneous  Techniques 

Multiple  Poisson  Distribution 

When  one  group  of  events  follows  the  Poisson  distribution  with 
the  parameter  mi  and  another  group  of  events  follows  the 
Poisson  distribution  with  the  parameter  m2,  the  population 
formed  by  combining  these  groups  will  follow  a  Poisson  dis- 
tribution with  the  parameter  (mi  +  ms).  (A  formal  mathe- 
matical proof  of  this  relationship  is  given  in  section  H  of  the 
Appendix.)  From  the  case  for  combining  2  groups  of  events  it 
follows  that  when  k  groups  of  events  are  combined  the  result- 
ing Poisson  distribution  has  the  parameter  (mi  -|-  m2  +  ms  + 
-f-  mfc).  One  traffic  application  of  this  relation- 
ship is  the  mathematical  description  of  the  total  arrivals  at  an 
intersection  when  the  arrivals  on  each  leg  are  known  to  follow 
the  Poisson  distribution.  Example  16  illustrates  this  application. 


so 
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Example  i6a 
Arrival  Rate— Southbound  (Durfee  Ave.)^ 


Number 

of  cars 

arriving 

Observed 

Theoretical 

per  ^o-sec. 

frequency 

frequency 

f 

interval 

f 

F 

F 

o 

i8 

15.0 

21. '6 

1 

32 

31.2 

32.8 

2 

28 

324 

24.2 

3 

20 

22.5 

17.8 

4 

13 

11.7 

14.4 

5 

1 

4-9 

6 

0 

1.6 

7 

1 

-9 

0.5 

-7.2 

11.3 

8 

1 

0.1 

^9 

0 

0.1 

120 


120.0 


122.1 


(Brackets  indicate  grouping  for  ^^  test.) 


120 


=  2.08 


=  0.125 


I 


P 


—  n  =  122.1  —  120.0  =  2.1 


u  =  6  —  2  =  4 


X%5  =  9-49 


16  See  footnote  9. 
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Example  16b 
Arrival  Rate— Northbound  (Durfee  Ave.) 


Number 

of  cars 

arriving 

Observed 

Theoretical 

per  ^O'Sec. 

frequency 

frequency 

f 

interval 

f 

F 

F 

0 

9 

5.6 

14-5 

1 

16 

17.2 

14-9 

2 

30 

26.3 

34-2 

3 

S2 

26.9 

18.0 

4 

19 

20.6 

17-5 

5 

10 

12.6 

7-9 

6 

S 

6-5l 

7 
8 

7 
3 

=14 

2.8 
1.1 

10.8 

18.1 

^9 

1 

0.4 

120  120.0  125.1 

(Brackets  indicate  grouping  for  the  x^  test.) 


368 
120 


3-07 


*'»  =  0.047 


I 


P 


■n  =  125. 1   —  120.0  =  5.1 


y  =  7  _  2  =  5 

X^05    =     I  I -I 
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Examp! 

e  16c 

Arrival  Rate 

—Northbound 

AND  Southbound  C 

ombine 

(Durfee 

Ave.) 

Number 

of  cars 

arriving 

Observed 

Theoretical 

per  ^o-sec. 

frequency 

frequency 

f 

interval 

f 

F 

F, 

o 

1 

o.f 

1 

4U4 

3-6M3-5 

14-5 

2 

9. 

9-2 

3 

19 

15-9 

22.7 

4 

20 

20.4 

19.6 

5 

16 

21.0 

12.2 

6 

19 

18.0 

20. 1 

7 

10 

13.2 

7.6 

8 

13 

8.5 

19-9 

9 

4l 

4-9 

lO 

1 

2-5 

11 

3 

9 

1.2 

9-5 

8.5 

^12 

1 

0.9 

120 


120.0 


125.1 


618 


=  5-15 


I 


f 


120 
"*  =  0.006 

■n  —  125. 1  —  120.0 
i;  =  8  —  2   =  6 
15-5 


rrin  +  ms  =  3.07  +  2.08  =  5.15 


5-1 


Xls 


Empirical  Formulas  Based  on  the  Poisson  Distribution 

Because  of  differences  between  theoretical  and  physical  con- 
ditions it  is  frequently  necessary  to  develop  empirical  formulas 
to  adequately  represent  the   true   behavior  of  traffic.   Where 
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traffic  arrives  in  a  random  manner,  the  Poisson  distribution 
may  serve  as  a  building  block  for  such  formulas. 

Consider  the  delays  at  urban  stop  signs  reported  by  Raff.^''^  If  a 
gap  of  t  seconds  is  required  for  crossing  the  main  street,  and  if 
the  side  street  cars  arrive  at  random,  the  percent  of  side  street 
cars  which  cross  without  delay  will  be  equal  to  the  percent  of 
gaps  in  the  main  street  traffic  which  are  t  or  greater.  From  the 
equation  B  this  is: 


100^ 


vt 

3600 


The  percent  of  side  street  cars  delayed  will,  therefore,  be 

IOOVI-^~3600J 

Raff  has  found  empirically  that,  because  of  sluggishness  of  the 
side  traffic,  the  delays  should  be  much  greater  than  the  above 
value  at  high  volumes  on  the  side  street.  He  has  developed  the 
following  empirical  formula  to  represent  this  situation. 

g -2.5Vg/3600  ^-2W/3600  "j 

I  00  <  I     — 


I   _  g-i/oV^/zma  Cj   _  g~vt/zwQ\ 


where  Vs  is  the  side  street  volume. 

When  Vs  =  o,  this  reduces  to  the  original  relationship. 

Combination  of  Poisson  Distribution  and  Other  Relationships 

Several  writers  have  used  the  Poisson  distribution  in  com- 
bination with  other  relationships  to  form  theoretical  expressions 
of  various  traffic  problems.  Although  these  techniques  are 
beyond  the  scope  of  the  present  discussion,  the  reader  may  wish 
to  familiarize  himself  with  some  of  the  papers  on  the  subject. 
The  following  partial  list  of  references  is  suggested  as  a  start- 
ing point  in  the  study  of  these  advanced  techniques. 

Garwood,  F.,  "An  Application  of  the  Theory  of  Probability  to  the 
Operation  of  Vehicular-Controlled  Traffic  Signals,"  Journal  of 

17  Raff,  Morton  S.,  A  Volume  Warrant  for  Urban  Stop  Signs,  The  Eno 
Foundation  for  Highway  TraflSc  Control,  1950. 
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The  Royal  Statistical  Society,  Supplement,  vol.  7,  no.  1,  1940, 
pp.  65-77. 

Raff,  Morton  S.,  "The  Distribution  of  Blocks  in  an  Uncongested 
Stream  of  Automobile  Traffic,"  Journal  of  The  American  Statis- 
tical Association,  vol.  46,  no.  253,  March  1951,  pp.  114-123. 

Kendall,  David  G.,  "Some  Problems  in  the  Theory  of  Queues," 
Journal  of  the  Royal  Statistical  Society,  Series  B,  vol,  13,  no.  2, 

195I'  PP-  151-173- 
Tanner,    J.    C,    "The   Delay    to    Pedestrians    Crossing   a    Road," 

Biometrika,  vol.  38,  parts  3  and  4,  Dec.  1951,  pp.  383-392. 

Pearcy,  T.,   "Delays  in  Landing  of  Air  Traffic,"  Journal  of  The 

Royal  Aeronautical  Society,  vol.  52,  Dec.  1948,  pp.  799-812. 

Statistics  of  the  Poisson  Distribution 

In  section  F  of  the  Appendix  the  mean,  variance,  and  stand- 
ard deviation  of  the  Poisson  distribution  are  derived.  The 
results  are: 

Mean:  jjlx  =  m 

Variance:  ax^  =  rn 

Standard  Deviation:  ax  =  y/m 

For  large  values  of  m  the  section  of  the  probability  curve  be- 
tween fij;  —  Ox  and  [Xx  +  Ox  is  approximately  68%  of  the  total 
area;  the  section  between  ix^  —  20-3.  and  ix.^  +  2ct^  is  approxi- 
mately 95%  of  the  total  area. 

Methods  of  Computation 

Review  of  Procedure 

The  following  is  a  review  of  the  procedure  of  computing 
probabilities  from  the  Poisson  distribution: 

I.  Determine  the  parameter  m.  This  parameter  is  the  average 
number  of  occurrences.  It  may  be  determined  from  observed 
or  assumed  data.  A  trial  may  consist  of  an  instantaneous  ob- 
servation, counting  events  during  a  time  interval,  counting 
events  in  a  unit  area,  etc. 

Total  number  of  events  observed 


m  = 


Total  number  of  trials  or  time  intervals,  etc. 
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2.  Once  m  has  been  determined,  the  probability  of  x  events 
occurring  at  any  trial  (during  any  time  interval)  is  computed 
from  the  formula 


W      P 

P{x)  =  "^- 


where  a:/  =  x(x—i)  (x  —  2)  —  3'2-i 

Table  i  lists  several  sources  of  tables  of  e"",  e~'^,  P{x),  etc. 
Figure  2  is  a  chart  of  the  cumulative  Poisson  distribution,  and 
may  be  used  for  rapid  calculations  where  accuracy  is  not  im- 
portant. A  method  of  making  slide  rule  calculations  is  indicated 
below. 

Chart  of  the  Cumulative  Poisson  Distribution 

Figure  2  is  a  modification  of  charts  by  Thorndike^^  and 
Working.i^  The  probability  P(x  ^  c)  is  plotted  against  m  with 
c  as  a  parameter.  The  probability  of  x  equal  to  or  less  than  c 
is  simply  read  from  the  curve.  To  obtain  P(x  =  c),  the  prob- 
ability of  exactly  c,  values  are  read  from  the  chart  for  x  =  c 
and  for  x  ^  (c  —1).  The  following  relationship  is  then  used: 

P(x  =  c)  =  F(x  ^  c)  —  P{x  ^  c  —  i) 

Slide  Rule  Calculations 

The  value  of  e"^  can  be  found  by  means  of  a  log-log-duplex 
slide  rule;  e~^  is,  of  course,  the  reciprocal  of  e'^.  Having  ob- 
tained e~'"  by  slide  rule  or  from  tables,  the  individual  terms  of 
the  Poisson  distribution  may  be  obtained  by  means  of  the 
following  relationship  for  P{x  -\-  1)  which  is  particularly 
adapted  to  slide  rule  calculations. 


P(x)  =  ^ 


18  Thorndike,  op.  cit. 

19  Working,  Holbrook,  A   Guide  to   Utilization  of  the  Binomial  and  Poisson 
Distributions  in  Industrial  Quality  Control,  Stanford  University  Press,  1943. 
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^x+lf,—m 


J-  \ 

■^  1 

^'        (x+i)l 

x/ix+i) 

;c+i 

Thus, 

it 

follows  that: 

P(o)  = 

^-m 

P(l)  = 

m 
I 

P(o) 

P(2)    = 

m 

2 

P(i) 

P(3)  = 

m 
3 

^2) 

P(x) 


etc. 

Table  i 
Aids  to  Computation  of  Poisson  Distribution 


Type  of  Aid 
and  Scope 

Source 

e^,  e-^ 
Tables   x   from   o.oo   to 
5.50  by  0.01  steps,  x  from 
5.0  to  10.0  by  0.1  steps 

Table  of  Exponentials 

Handbook     of     Chemistry     and 

Physics 

Chemical       Rubber      Publishing 

Company,  Cleveland 

Various  editions 

,      Tables 
xl 

Individual    and    Cumu- 
lative 

a  from  0.00 1  to  100.0  by 
various  steps 

Molina,  E.  C,  "Poissons  Exponen- 
tial  Binomial   Limit,"   Van   Nos- 
trand,  1945 

— : —  Tables 

V- 

e  from  0.1  to  1.0  by  steps 
of  0.1.  €  from  1.0  to  20.0 
by  steps  of  1.0 
Individual  and  cumula- 
tive 

Appendices  VI  and  VII 

Fry,   T.  C,   "Probability  and  Its 

Engineering  Uses,  Van  Nostrand, 

1928 
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Type  of  Aid 
and  Scope 

Source 

m'^  e-'^ 

Appendix  VI 

Greenshields,   Bruce  D.;   Shapiro, 
Donald;      Ericksen,      Ehoy      L.; 
"Traffic    Performance    at    Urban 
Street    Intersections,"     Technical 
Report    No.    1,    Yale    Bureau    of 
Highway  Traffic,  1947. 

x! 
Cumulative    Tables     m 
from  0.1  to  15.0  by  steps 
of  0.1 

Chart 
Probability  that  no.  cars 
arriving  in  a  given  time 
will   be   X   or  more   (or 
less  than  x) 

Figure  3 
Adams,  William  F.,  "Road  Traffic 
Considered  as  a  Random  Series," 
Institution     of     Civil     Engineers 
Journal,  Nov.,  1936,  pp.  121-130-I-. 

Chart 

Cumulative 

P(x  ^  c) 

Figure  5 
Thorndike,  Frances,  "Applications 
of    Poissons    Probability    Summa- 
tion,"    Bell     System     Technical 
Journal,  vol.  5,  No.  4,  Oct.,  1926, 
pp.  604-624. 

Slide  Rule 

X  from  0.001  to  10.0  on 
3  or  4  10-inch  scales  plus 
normal    slide    rule    fea- 
tures 

The  Fredrick  Post  Company 

Versalog  Slide  Rule 

X  from  0.001   to  10.0  on  4  scales. 

Keuffel  and  Esser  Company 

Log-Log-Duplex  Trig  Slide  Rule 

Log-Log-Duplex     Decitrig     Slide 

Rule 

X  from  0.01   to  10.0  on  3  scales. 

e^  and  e~* 
Tables  x  from  0.0  to  0.1 
by  steps  of  0.0001;  from 
0.1    to    3.0    by   steps   of 
0.001;  from  3.0  to  6.3  by 
steps  of  0.01;  from  6.3  to 
10  by  steps  of  0.1 

Hayashi,     Keiichi,     "Fiinfstellige 
Tafeln  der  Kreis-  und  Hyperbel- 
funktionen  sowie  der  Funktionen 
e^  und  e-^"  Berlin,  1944,  Walter 
de  Gruyter  &  Co. 
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Gaps  oil  Arroyo  Seco  Freeway 

Lane  1,  2  to  2:30  P.M. 

October,  1950 
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Greek  Alphabet 

Name 

Lower  Case 

Upper  Case 

Alpha 

a 

A 

Beta 

13 

B 

Gamma 

y 

r 

Delta 

8 

A 

Epsilon 

€ 

E 

Zeta 

C 

z 

Eta 

V 

H 

Theta 

e 

© 

Iota 

I 

I 

Kappa 

K 

K 

Lambda 

A 

A 

Mu 

/* 

M 

Nu 

V 

N 

Xi       ^' 

$ 

S 

Omicron 

0 

o 

Pi 

TT 

n 

Rho 

P 

p 

Sigma 

U 

S      - 

Tau 

T 

T 

Upsilon 

V 

Y 

Phi 

<i> 

* 

Chi 

X 

X 

Psi 

^ 

^ 

Omega 

0) 

o 
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Fundamental  Concepts 

Permutations  and  Combinations 

The  subject  of  permutations  may  be  easily  illustrated  by  the 
example  of  code  words.  Consider  that  3  cards  marked  A,  B,  C 
are  available  and  are  to  be  used  to  form  as  many  3-letter  code 
words  as  possible.  The  result  will  be: 

/B  C  ABC 

K 

\C  B  A  C  B 

/AC  B  A  C 

b/ 

\C  A  B  C  A 

,A  B  CAB 

\b  A  C  B  A 

In  each  case  there  will  be  3  choices  for  the  first  letter,  2 
choices  for  the  second  letter,  and  1  choice  for  the  third  letter. 
From  this  it  follows  that: 

Pn  =  Permutations  of  n  things  taken  n  at  a  time 

=  n{n—i)  (n  —  2)  •  •  •  •  3  •  2  •  I 

=  n/  =  factorial  n 

Pf  =  3  •  2  •  I   =6  which  checks  the  empirical  result 

If  there  are  5  cards  marked  A,  B,  C,  D,  E  and  it  is  desired 

■  that  3-letter  code  words  be  formed,  there  will  be  5  choices  for 

the  first  letter,  4  choices  for  the  second,  and  3  for  the  third: 

P3  =  Permutations  of  5  things  taken  3  at  a  time 

=  5  .  4  .  3  =  60 

_  5  •  4  •  3  •  2  •  I  _  5:  _       5! 

2-1  2!  (5-3)! 
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For  the  general  case  of  the  permutations  of  n  things  taken  m 
at  a  time 


{n  —  m)i 


In  these  code  words  the  order  is  important,  for  ABC  is  different 
from  ACB.  Consider,  however,  3-card  hands  formed  from  the  5 
cards  A,  K,  Q,  J,  10.  If  order  were  important,  there  would  be 
60  permutations  as  follows: 

AKQ   AQJ  AJio    KQJ  AKJ  AKio  AQio    KJio  KQio  QJio 

AQK    AJQ   AioJ    KJQ   AJK  AioK  AioQ   KioJ  KioQ    QioJ 

KAQ   QAJ   JAio    QKJ    KAJ  KAio  Q.A10   JKio  QKio    JQio 

KQA    QJA   JioA  QJK    KJA  KioA  Q.10A   JioK  QioK   JioQ 

QAK   JAQ    loAJ   JKQ  JAK  loAK  loAQ    loKJ  loKQ,   loQJ 

Q.KA   JQA    loJA   JQK  JKA  loKA  loQA    loJK  loQK    loJQ 

Note  that  each  vertical  group  is  P|  while  the  whole  array  is 
Pi.  In  general,  however,  the  order  of  the  cards  is  unimportant. 
It  is  the  particular  group  of  cards  that  is  important  (represented 
by  the  10  groups  above).  When  order  is  unimportant  the  group- 
iuCT  is  known  as  a  combination. 


^5  ^  Combinations  of  5  things   ^  Permut.  of  whole  array 
^  taken  3  at  a  time  Permut.  within  combination 

P|        60 
=  Pi  =  "6-  =  ^° 


pre 


PZ        (n-m)/  m! 
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Factorial  Zero 

For  clarity  in  certain  problems  it  has  been  found  convenient 
to  define  the  factorial  zero  as  follows: 

n' 

(n-i)  !  =- 

n 

let:  n  =  I 

I ' 
then:  (i  — i)  !  =  — 

o  !  =  I 

Laws  of  Probability 

The   following  are   two   important   laws   concerning  prob- 
ability: 

1.  Total  Probability 

If  two  events,  A  and  B,  are  mutually  exclusive  (if  A  occurs, 
B  cannot  occur  and  vice  versa)  the  total  probability  that  one 
of  these  events  will  occur  is: 

P(A  or  B)  =  F(A  +  B)  =  P(A)  +  P(B) 

2.  Joint  Probability 

If  two  events,  A  and  B,  are  independent  (the  occurrence  of 
one  has  no  influence  on  the  other)  the  probability  that  both 
will  occur  together  is: 

P(AB)  =  P(A  and  B)  =  P(A)P(B) 


APPENDIX  C 


The  Binomial  Distribution 

Consider  a  population  in  which  each  item  may  possess  one  of 
two  mutually  exclusive  characteristics  (head  or  tail,  good  or 
bad,  o  or  1,  etc.) 

Let:  p  =  probability  of  occurrence  of  characteristic  A 

^  ^  (i  —  P)  =  probability  of  occurrence  of  characteristic 
B  (non-occurrence  of  A) 
Then,  by  the  law  of  total  probability, 

P(A  or  B)  =  p  -f  ^ 

Suppose  that  a  sample  of  n  items  is  drawn  from  the  popula- 
tion under  the  following  conditions: 

a.  The  size  of  the  population  is  infinite.  (This  restriction  insures 
that  withdrawing  the  sample  does  not  alter  the  relative  pro- 
portion of  A  and  B  remaining  in  the  population.  The  same 
result  may  be  achieved  with  a  finite  population  by  drawing 
one  item,  replacing,  stirring,  drawing  the  next  item,  etc.) 

b.  The  sample  is  selected  from  the  population  at  random. 

Then,  by  the  law  of  joint  probability: 

the  probability  that  all  n  items  in  the  sample  are  A's  =-p'* 


the  probability  of  (n — i)  A's  and  i  B  =  p'^-'^q 
the  probability  of  {n — 2)  A's  and  2  B's  =  p^-^q^ 

etc. 
the  probability  of  lA  and  (n — 1)  B's  =  pq^-'^ 
the  probability  of  n  B's  =  q"' 


r     i 


where  the 
order  of 
A's  and 
B's  is  sig- 
nificant 


In  general,  for  m  A's  and  (n—m)  B's  the  probability  =  p"^  q"--"* 
where  m  ^  o,  1,  2,  •  .  .  (n — 1),  n 

In  drawing  the  sample,  the  order  in  which  the  A's  occur  can 
take  on  many  possibilities,  the  number  being  equal  to  the  com- 
binations of  n  things  taken  m  at  a  time  =  C^ 

45 


46  POISSON    DISTRIBUTION 

For  instance,  if  the  number  of  items  in  the  sample  is  5,  there 
are  Cf  =  10  ways  in  which  a  sample  composed  of  gA's  and  2B's 
may  be  drawn.  Each  of  these  ways  will  have  a  probability  of 
p^q^.  Thus,  by  the  law  of  total  probability,  the  probability  of 
3  A's  and  2  B's  is: 

I  op^q^ 
or  in  general: 

C^pmqn-m 

Considering  the  various  possibilities  of  o,  1,2,  .  .  .  .n  A's  and  n, 
(n— 1)  ....  1,  o  B's,  the  total  probability  of  occurrence  of  A's 
and  B's  is  given  by 

n 
P{A,  B)  =      I,     Clp'^q''-^ 
m  =  o 

The  right  hand  member  is  equal  to  (q  +  p)"  and  hence  this 
relationship  is  known  as  the  binomial  distribution. 

n 

m  =  o 
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Derivation  of  the  Poisson  Distribution 

Introduction:  The  Poisson  Distribution  is  applicable  to  popu- 
lations having  the  following  properties: 

a.  The  probability  of  occurrence  of  individuals  having  a  par- 
ticular characteristic  is  low. 

b.  The  characteristic  is  a  discrete  variable. 

The  Poisson  Distribution  as  a  Limiting  Case  of  the  Binomial 
Distribution 

Let    n  =  number  of  items  in  sample 

p  =  probability  of  occurrence  of  a  particular  characteristic  E 

^  =  (i  —p)   —  probability  of  non-occurrence  of  character- 
istic E 

X  =  number  of  items  in  sample  having  characteristic  E. 
Then,  from  the  binomial  distribution: 

P(x)  =  C^p^q""-^  =  c^p^ii-py-' 

X  =  o,  I,  2  •  •  '  •  n 
Now  let: 

p  be  made  indefinitely  small 

n  be  very  large 

pn  =  m,  where  m  is  finite  and  not  necessarily  small. 


Then:  p  =  — 
n 
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_  n!        (  ^\  (     _    ^  V~^ 

x/(n  —  x)/  \  «  /     \  "  / 

_  nf        (  ^\  (  '^ \"  /  _    ^ ^ 

x!{n  —  x)!  \n  )    \  n  J    \  n  J 

(n  —  x):  rf'  \  1 | 

L  \       «/  J 


=  U]  [5]  [C] 


Now,  if  w— >  tx> 
/iw  P{x)  =  Urn 


JM[5][C] 

_  r/zm  A  1  riim  B  1  f/zV/i  C  1 

1^;?— »ooJ    |_/2— ^ooj    |_?z— »c»J 


A    = 


n— >  00  x/ 


-  =  (■-")" 


/zm  B    =  e-^ 


(See  appendix  E  for  proof) 


n— >oo 


C  = 


(n  — x) 


■'"'  (■  -  ^)" 
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When  n  is  verjr  large,  negligible  error  is  introduced  by  repre- 
senting n!  by  one  term  of  Stirling's  formula.  The  same  state- 
ment holds  for  (n—x)\ 


.-.  C= 


c  = 


V2: 


irn     n"     e 


n       a—n 


■\/  27r  {n  —  x)  (n  —  xy-''  ^-("-^)  (i   —  — )    "" 

■\/27r(ra  — a:)^  (n  —  x)"-''  ^-("-^)  (i   —  —  )    «'' 

(n  — A 


("- 


-)"-'  ('  -  f)' 


[-]  [■  ^  ;I 


(■  -  r  (■  ^  ^)- 


H  [■  -  r 


(■  -  0-  _(■  -  ^) 


[c]  [c]  [o]  [c] 


Urn  Ci  =  ^~* 


(■-3 


a  =   I  -  - 


r-i 
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lim  Ci  =  I 
n— ^00 

C,  = 


{■  ^  0' 


lim  C3  = 
n—>co 


(See  Appendix  E) 


C,= 


/  mY 


lim  C4  =  I 

n— >oo 

l^^  Q    _  I  /z>?2  Cil  Vlim  Ci]  Vlim  C3 
^~^  °°         I  n— >  00  J  I  n— >  00  J  I  n^  00 

-  [-]  [  ■  ]  M" 

=   I 

km  P{x)  =  —^  e-^ 


lim  C4 
n — >oo 

I 


Since  the  main  body  of  this  discussion  assumes  the  existence  of 
the  conditions  for  the  Poisson  distribution  the  above  equation 
may  be  written  simply: 


^w  =  ^ 


APPENDIX  E 


('  -  ^y  s 


Derivation  of  Limit  I  i  ~~  —  I     (Used  in  the  derivation  of  the  Poisson 


n— >oo    \  ••- /     distribution) 


Let     n  =  - 

X 


Urn 
n-^co 


(■-^«)"-; 


lim    (i  —mx)'^;  where  n  = 


=  lim 


I 


(t) 


(mx)  + 


(t)G  - '}'"''' 


2! 


>  00 

X 


m  -  %  -  ^y-'- 



3!      ^ 

lim      \  —  m -\ -Ai—x) -^{i—x){i—ix)-\- 

I 

—  — ^00 



x 

,    n?       m^    , 
=  I  -^  +  ^-^  + 

But  expanding  e~^  in  a  McLaurin  series  gives: 


,    III         III,     , 
e-^  =    I    -  m  +   -r ;-  + 


2!    31 


lim 


{-")- 
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Mean,  Variance,  and  Standard  Deviation  o£ 
Poisson  Distribution 

Mean 

The  mean  of  any  continuous  function  is  obtained  by: 

^   /  xf(x)dx 


f  f{x)dx 

(The  mean  may  be  considered  as  the  distance  to  the  center  of 

gravity.) 

For  discrete  functions  the  comparable  function  defining  the 

mean  is: 

^  ^xf{x) 
When  dealing  with  probabilities: 

m  =  PM 

ff{x)dx  =   I 


S     P(x)  =  I 
x  =  o 


ix  =   fxP(x) 


fx  =     "Ex  P(x) 
x  =  o 
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F 

For  the  Poisson  distribution 

P(x) 

M 

.^0      -' 

=  o  +  me^""  + 

2! 

-  + 

3! 

=  me~^     I  +  m 

^  2! 

=  me~^e"^ 
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By  definition,  the  variance  a^  may  be  expressed: 

2/W  (^-m)- 


c"  = 


2/W 
For  the  Poisson  distribution 

'E(x-my  P(x)        'E{x^-2xm-{-m'^)  P(x) 


a    = 


i:p(x)  I 

=  'Ex^Pix)  -  2mi:xP(x)  +  m^i:Pix) 
The  last  two  terms  reduce  as  follows: 

—  amExPix)  =  —2m(m)  =  —2rr? 
m^1.P{x)  =  m2 
The  first  term  may  be  reduced  by  the  following  steps: 
i:x^P(x)  =  2  [x{x-i)  +  x]  P(x)  =  SAr(x-i)  P(x)  +  i:xP{x) 
=  Sx(x-i)  P(x)  +  m 


=  To  +  o 


2!  q!  4-! 
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-—-      m'^    D     ^ 


(i+ZTzH -y  -\-  •  •  •   \  -\-  m 


cr^  =  m?  -{-  m  —  irr?  -\-  m^  =  m 

Standard  Deviation 

By  definition  the  standard  deviation,  a,  is  expressed  by 

a  =  V^ 
Tlius,  for  the  Poisson  distribution 

a  =  V  m 


APPENDIX  G 

Notes  on  the  -^  Test  for  Goodness  of  Fit 
Let:   /  =  Observed  frequency  for  any  group  or  interval 

F  =  Computed  or  theoretical  frequency  for  same  group 
Then,  by  definition: 

x^^    s    ili^-Il.  (I) 

where  k  =  number  of  groups 
Expanding: 

^     f^        if-  F         F^ 
■         F-  F  F- 

k       r  2  k  k 

=     S    -^  -  2     S  /i  +     S   Fi 

1=1     *  i=i  i=i 

But  by  assumption  in  the  fitting  process: 

k  k 

S    fi=     S     Fi  =  n 
1=1  i=i 

where  n  =  total  number  of  observation 
So  that: 

Either  equation  (i)  or  equation  (2)  may  be  used  for  purposes 
of  computation.  Usually  (2)  will  simplify  the  amount  of  work 
involved. 
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The  value  of  x^  obtained  as  above  is  then  compared  with  the 
value  from  Figure  3  or  from  tables  of  x^-  Such  tables  may  be 
found  in  any  collection  of  statistical  tables,  and  relate  the 
values  of  x^  and  significance  level  with  the  degrees  of  freedom. 
The  number  of  degrees  of  freedom,  v,  may  be  expressed:^"- ^^ 

V  =  G— 2 
where  G  =  number  of  groups 

For  this  value  of  v  to  be  valid,  however,  it  is  necessary  that  the 
theoretical  number  of  occurrences  in  any  group  be  at  least  5. 
One  writer^^  further  stipulates  that  the  total  number  of  observa- 
tions be  at  least  50.  When  the  number  of  theoretical  occur- 
rences in  any  group  is  less  than  5,  the  group  interval  should  be 
increased.  For  the  lowest  and  highest  groups  this  may  be 
accomplished  by  making  these  groups  "all  less  than"  and  "all 
greater  than"  respectively. 

20  Dixon,  W.  J.,  and   Massey,  F.  J.,  Jr.,  Introduction   to  Statistical  Analysts, 
McGraw-Hill,  1951,  pp.  190-191. 

21  Yule,   G.   Udny,   and   Kendall,   M.   G.,  An  Introduction   to   the   Theory  of 
Statistics,  14th  edition,  Haffner  Publishing  Company,  1950,  p.  475,  p.  476. 

22  Cramer,  Harald,  Mathematical  Methods  of  Statistics,  Princeton  University 
Press,  1946,  p.  435. 
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Derivation  of  the  Distribution  of  the  Sum  of  Independent 
Poisson  Distributions 

Consider  a  population  made  up  of  two  subpopulations  A  and  B, 
each  distributed  according  to  the  Poisson  distribution. 
For  subpopulation  A 

X       — m 
-y,  a  „       a 

P{Xa)    =    "" 


For  subpopulation  B 


h  „       b 
PM     =    ■■■ 


— m, 

m  "  e 


Xb! 


If  k  items  occur  in  a  trial  from  the  total  population,  there  may 
be  a  mixture  of  Xa  and  x?,  as  follows: 


Xa    —  kj 

^6  = 

o 

Xa   -\-   Xb    =    k 

Xa    =    k    — 

I 

Xb    = 

I 

Xa   -\-   Xb    =    k 

Xa    =    k    — 

2 

Xb    = 

2 

Xa   -\-   Xb    =    k 

Xa    =    2 

Xb    = 

k  - 

2 

Xa   +   Xb    =    k 

Xa    =    I 

Xb    = 

k  - 

I 

Xa   -\-   Xb    =    k 

Xa    =    O 

Xb    = 

k 

Xa   -\-   Xb    =    k 

P{k)    =    P{Xa    =    k,  Xb    =    O)    -\-   P{Xa    =    k    —    I ,  Xb    =    l) 

+   PiXa    =    I,  Xb    =    k    ~    l)    -\-   P(Xa    =    O,  Xb    =    k) 

— m       — m  — m  — m  — m  — m 

ma"  e     "■  e     °        rria      e     "■  rrib  e     "        m/     e     °-  m^  e     ° 


kl  0\  '  (k-l)\l\  '  (k-2)\2\ 
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— m  — m  - 

nia  e     "  mfe"""^  e     °    ,    e     "■  rrib   e 


a  *,o,fc— 1   „       b 


_^   ___,^^,__    _|_ 


i!(^-i)!  '         o\  k/ 


p/,x   _    -^<j   -^j  jmj'        k  mj'  ^  nib    ,    k{k—i)  iria^  ^  /nb^ 
^^^^  ~  '        '       \k!    "^    k{k-i)\    "^   yt(A;-i)  (^—2)!  2 

A:(A;—  i)  •  •  •  3  ■  2  •  I  ?/2a  r^ft''"^    ,    rnb'\^ 
k{k-i)  •  •  •  3  •  2  •  I  (/t-i)!    ^17/ 


+ 


P{k)  =  — <ma''  +  krria''  ^  rtih  i — 

+  •  ■  ■  •  -\-  k  rria  mb^~'^  +  rub^ 

—  (m  +  TO  ) 


} 


P(^)  =  : ^J^^^  +  mf 

When  there  are  subpopulations  A,  B,  .  .  .  .  ,  ■Z,hy  application 
of  a  similar  argument  the  distribution  for  the  whole  population 
is  found  to  be 

1  I  M.       "("*       +    TO       +••••+    TO    ) 

p,j.   ^   (ffig  +  ma  +  •  •  •  •  +  m^Y  e      °        " ^_ 

k! 
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A  theoretical  study  of  conditions  affecting  the  traffic  of  vehicles 
on  a  highway  requires  the  constant  use  of  probability  theory. 
On  the  one  hand,  the  distribution  of  vehicles  in  each  lane  is 
in  part  a  matter  of  chance.  On  the  other  hand,  whenever  one 
studies  the  behavior  of  any  large  number  of  individuals,  in- 
evitable departures  from  the  laws  which  apply  to  their  totality 
are  found  by  statistical  analysis  to  follow  certain  empirical  and 
relatively  stable  patterns.  These  are,  in  effect,  special  laws  of 
probability. 

Reprinted  with  permission  from  Travaux,  January  1955,  revised  and  updated 
by  the  author.  Translated  from  the  French. 
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In  particular,  Poisson's  law,  applying  to  rare  events,  has  thus 
far  been  the  chief  theoretical  instrument  for  dealing  with  prob- 
lems of  vehicular  traffic  on  two  or  three  lane  highways,  espe- 
cially the  problem  of  determining  the  distribution  of  such 
traffic  both  in  time  and  in  space.  This  law,  as  is  well  known, 

assigns  the  value  e  ~^^  - — p-  to  the  probability  that  there  are 

"n"  vehicles  in  a  time  interval  0  chosen  at  random,  N  being 
the  average  number  of  vehicles  in  unit  time.  The  probability 
that  there  are  no  vehicles  in  this  interval  is 

g-NB 

While  the  results  yielded  by  this  formula  agree  well  enough 
with  actual  observation  when  the  traffic  density  is  low  (a  few 
dozen  cars  per  hour  per  lane),  they  differ  widely  from  reality 
when  the  density  is  significantly  larger.  The  reasons  for  this 
discrepancy  between  theoretical  prediction  and  the  observed 
data  are  clear  enough.  But  without  going  into  any  analysis  of 
physical  causes,  it  seems  obvious  that  the  gap  between  theory 
and  actuality  can  be  considerably  narrowed.  It  is  the  object  of 
this  paper  to  show  how  this  can  be  done. 

We  shall  suppose  that  the  entire  set  of  spacings  between  suc- 
cessive vehicles  consists  of  a  number  of  distinct  parts  or  sub- 
sets, each  having  distinct  mean  values  and  each  obeying  some 
Poisson-type  law.  For  simplicity  we  consider  just  two  sub-sets. 
Let  the  number  of  spacings  per  unit  time  for  these  two  sets  be 
respectively 

YJVand  (i  —  7)JV 

with  mean  spacing-values  of  ti  and  t2  seconds  respectively,  and 

with  ti  <  ^2-  We  also  suppose  that  the  entire  set  of  spacings  is 

a  set  of  random  and  independent  elements. 

Then  the  probability  that  there  is  no  vehicle  in  an  interval 

_±  _± 

6  would  be  e     *i  if  all  spacings  were  in  the  first  sub-set  and  e     '2 

if  they  were  all  in  the  second  sub-set. 
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Now  the  respective  times  covered  by  the  two  sub-sets  are 

7 jWi  and  i  —  yNh, 
and  clearly 

i—yNh  =  {i—y)J^k 

Accordingly,  the  probability  that  there  are  no  vehicles  in  an 
interval  6  will  be  given  by 

Pie)  =  Nyhe     "   +  ii-y)Nt2e     ''  (i) 

with  the  relation 

Nyh  +  ^1-7)^2  =  I, 

which  merely  states  that  the  N  spacings  cover  unit  time.  Before 
proceeding  further,  it  must  be  observed  that  the  first  set  of 
spacings  might  apply  to  retarded  vehicles  which  are  prevented 
from  passing  by  opposing  traffic,  and  the  second  set  to  free- 
moving  vehicles  which  are  able  to  pass  at  will.  As  vehicles 
cannot  be  considered  as  mere  points,  two  successive  vehicles  in 
the  first  set  must  necessarily  be  separated  by  a  time  interval 
having  a  positive  lower  bound  e.  On  the  contrary,  free-flowing 
vehicles  having  opportunities  to  pass  may  exhibit  spacings 
equal  to  zero. 

Hence  the  law  of  spacings  which  we  use  in  practice  is  not 
formula  (i),  but  rather  the  formula 


P{d)  =  Ny{h-e)e     ''-'  +  N{i-y)he     ''  (2) 

We  must  now  examine  the  practical  usefulness  of  the  laws 
embodied  in  formulas  (1)  and  (2)  in  providing  answers  to 
various  questions  arising  out  of  the  study  of  traffic  problems. 
They  are  considered  in  some  detail  in  what  follows. 

First  Question:  Probability  of  a  spacing  longer  than  x. 

Suppose  that  the  vehicles  in  our  set  are  represented  by  points 
on  a  straight  line  and  that  p{x)dx  is  the  probability  of  a  spac- 
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ing  between  x  and  x  -\-  dx.  What  then  is  the  probability  P{9), 
that  an  interval  6  chosen  at  random  on  the  straight  line,  con- 
tains no  points,  all  choices  of  6  on  the  line  being  equally 
probable? 

Now  the  probability  that  the  initial  point  A  of  the  interval 
AB  =  6  will  be  in  a  spacing  lying  between  x  and  x  -{-  dx  is 

kxp(x)dx, 
k  being  defined  by  the  condition 


/ 


kxp{x)dx  =  I; 


which  simply  states  that  the  point  A  necessarily  falls  in  a  spac- 
ing which  lies  between  O  and  oo.  The  value  of  k  is  given  by 


/ 


^■, 


xp{x)dx 


where  Xm  is  the  mean  value  of  x. 

If  the  interval  AB  =  0  contains  no  points,  L  being  the  first 
point  preceding  A,  and  M  the  first  one  following  A,  with 
LM  =  X,  it  follows  that  x  ^  9  and  that  the  initial  point  A  lies 
between  L  and  P,  where  LP  =  x  —9  ^  O.  Now  the  probability 
that  the  interval  AB  =  9  contains  no  points  is  the  same  as  the 
probability  that  its  initial  point  A  lies  between  L  and  P^  and 
this  probability  is  found  from  the  above  to  be 


('-?)' 


:p(x)dx  =  N{x-9)p{x)dx 


^r- 


X- 


P 


B 


M 


Figure  i 
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Hence  upon  integration  over  all  values  of  x  between  9  and  co, 
we  have 

N{x-d)p{x)dx  =     I     N  dx    I    p{x)dx  (3) 

From  this  we  find  by  two  successive  differentiations,  that 

I     d'  P(x) 

Using  this  result,  we  find  by  taking  the  value  of  P{x)  from  (1) 
that  in  this  case 


p{x)  =  -j-e         +  -— ^  e  (4) 


y  <i     ,     I  —  7         h 

—  e         -\ e 

h  h 

while  the  total  probability  of  a  spacing  longer  than  x  is  given  by 

I  X  _x_ 

p(x)dx  =  ye~'^  +(i-y)e     ''  (5) 


/ 


If  on  the  other  hand  we  take  the  value  of  P(x)  furnished  by  (2), 
then 


p(x)  =  ~-e  +  -— ^  e  (6) 

'1  —  C  l2 


and 


=/ 


/)(x)6fA:  =  7^       '    '    +(1—7)^       '  (7) 


It  will  be  seen  by  examining  the  graphs  of  Figure  2  which  has 
been  taken  from  Statistics  with  Applications  to  Highway  Traffic 
Analyses  by  Greenshields  and  Weida,  that  there  is  good  agree- 
ment between  the  observed  data  and  those  furnished  by  formula 
(7).  The  legend  on  the  diagram  gives  the  actual  numerical 
values  of  the  various  parameters  used  in  this  particular  case. 
(See  Appendix,  p.  73.) 
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Figure  2 
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Figure  3  gives  another  example,  studied  on  a  French  road, 
showing  how  actual  data  are  approximated  by  the  theoretical 
formula.  There  were  312  spacings  observed  in  a  time  period  of 
52  minutes,  45  seconds.  Opposing  traflfic  amounted  to  only  12 
percent  of  the  total  flow. 
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w,     _ 
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—  Theoretical  resuhs 
n  =  n   +  n 


The  dotted  cun'es  give  the  distribution  of 
spacings  of  the  first  class  (constrained  vehicles) 
and  of  the  second  class  (free  moving  vehicles) 


n:  Number  of  spacings 

(cumulated) 
6:  Length  of  spacings 
in  seconds 


4  5 

Figure  3 


6" 


Second  Question:  Probability  that  an  interval  6,  taken  at 
random,  contains  no  vehicles,  but  is  bounded  by  a  vehicle 
on  one  side. 

First,  we  clarify  the  significance  of  the   iterated  integrals 
occurring  in  formula  (3). 
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If  AB  =  6  is  placed  at  random  on  the  spacing  LM 
can  be  empty  only  if  A  and  B  are  between  L  and  M. 


X,  AB 


Figure  4 

Accordingly,  x  ^  6,  and  Figure  4  shows  that  the  various  cases 
in  which  AB  contains  no  points,  may  be  obtained  by  letting  L 
range  over  all  positions  preceding  A,  and  M  over  all  positions 
following  B. 

If  B  lies  between  M  and  M',  where  MM'  =  dv,  and  if  kdv 
is  the  probability  of  this  happening,  then  the  probability  that 
LM  lies  between  u  and  u  -j-  du  is  the  probability  of  a  spacing 

u-{-d-\-v<x^u-\-du-\-d-{-v, 

that  is, 

p(u  +  d  -\-  v)du 

Therefore  the  total  probability  that  L  precedes  A,  M  being 
given  on  dv,  is 


k  dv 


+  0  +  v)du 


Also  the  total  probability  that  L  precedes  A[u'^  o)  and  that  M 
follows  Biv"^  o)  is 


f      kdv    I 

0  c/  o 


p{u  +  d  +  v)du 


If  we  make  the  substitution 

u  +  d  +  V 
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this  can  be  written  in  the  form 


/       k  dv     I        p(x)dx 

Jo  J  B+v 


By  the  further  substitution 

e  +  V  =  V, 

we  find 


k  dy     I     p{x)dx  =  k     j      p{x)dx    /       dy, 
J  y  Je  Je 


or 


/oo 


-6)p(x)dx 


where  k  is  determined  by  the  condition  that  P{0)  =  1,  that  is, 
k  =  N,  and  we  recover  formula  (3). 

The  upshot  of  this  argument  is  the  demonstration  that  there 
are  several  distinct  definitions  of  the  probability  that  a  given 
time  interval  is  empty.  In  particular  our  second  question  can 
now  be  answered,  for  the  probability  that  a  time  interval  be- 
tween 6  and  6  -}-  d6  is  empty  but  bounded  on  one  side  by  a 
vehicle,  is  found  from  the  above  results  by  taking  MB  =  v  =  0, 
while  L  of  course  precedes  A .  In  fact,  if  t/  =  O  we  obtain  from 
formula  (3)  by  differentiation 


-"^-Ldd  =  Ndd 

da 


/CO 
p{x)dx 


We  shall  hereafter  write 


Third    Question:   Probability    that    an    interval   6,    taken    at 
random,  contains  exactly  n  vehicles. 

If  Xi,   X2,  Xn  are  the  distances  of  the  various  points 
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(vehicles)  interior  to  6  measured  from  the  initial  point  A  of  the 
interval,  we  can  use  the  preceding  results  to  find  the  prob- 
abilities that  the  intervals  Axi,  X1X2, Xn-i  Xn  are  empty  but 

bounded  on  the  right  by  a  point  (vehicle),  namely  Xi  for  the 
first  interval,  X2  for  the  second,  and  Xn  for  the  last. 

Consider  the  interval  Axi.  The  probability  that  it  has  the 
required  property  is 


/ 


J{xi)dxi 


Now  the  point  Xi  must  be  followed  by  a  spacing  smaller  than 
6  —Xi,  that  is,  the  corresponding  interval  must  be  empty  but 
bounded  on  the  right  by  a  point,  viz.  X2.  The  probability  of  this 
happening  is 


I      Jixi)dxi     I        J(x2  —  Xi)dX2 


We  continue  in  this  way,  up  to  the  last  interval  B  —Xn,  which 
is  empty  but  not  bounded  on  the  right.  Hence  the  desired 
probability  is 

I       J(xi)dXi      I       J(x2  —  Xi)dX2 /       J(Xn  —  Xn-l)dXn  P(d  —  Xn) 

Jo  %J  xi  tJ  Xn—i 

(JV6)" 
It  is  this  formula  which  ffives  the  value  e  ~^^  ^ — /~  when 

n! 

P(e)  =  e  -^'. 

Fourth  Question:  Probable  delay  in  waiting  for  an  empty 
interval  9. 

This  question  comes  up  whenever  a  vehicle  in  one  lane  wishes 
to  pass  another  car  ahead  of  it,  also  whenever  a  vehicle  enters 
upon  or  cuts  across  a  main  highway.  We  first  determine  the 
probability  that  an  interval  of  length  I,  whose  initial  point  is 
chosen  at  random,  contains  no  void  of  length  6. 

This  is  obviously  a  function  F{1,  6)  which  steadily  decreases 
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as  I  increases,  0  being  given.  The  decrease  in  F{1,  6)  due  to 
an  increase  dl  is  the  probability  that  the  interval  /  -\-  dl  con- 
tains a  void  6,  while  /  itself  contains  no  void.  It  is,  in  other 
words,  the  probability  that  the  interval  /  contains  a  void  6  — dl 
at  its  extremity,  dl  itself  being  empty.  Consequently  it  is  the 
product  of  the  probabilities  that  /  +  dl  —6  contains  no  void  of 
length  6  and  that  9  —dl  is  empty  but  bounded  on  one  side  by 
a  vehicle,  provided  we  neglect  all  infinitesimals  of  order  higher 
than  that  of  F{1  —9,  e)j{9)dl.  We  can  therefore  write 

bF{l,  9) 


81 


=  -  F(l-9,  9)J9 


This  relation  enables  us  to  define  F{1,  9)  step  by  step  once  we 
know  it  on  an  interval  of  minimum  length  9. 

Now 

F{1,  9)  =  I         for         0  ^  I  <  9 
and  from  this  we  obtain 

8F 

^  =  -  J{9)    and   F  =  K-J{9)1  for    0  ^  I  <  29 

Since  obviously 

F{9,9)  =  i-P(9), 
we  see  that 

K=  i-P(9)  +J{9y9, 
and 

F{1,9)  =  i-P{9)  -J(9)(l-9) 
Continuing  in  this  manner  to 

n9  ^  I  <  (n-\-i)9 
we  find  that 

F(l,  9)  =  1  +  S(-i)"    ^^^  (l-n9)-^  +  ^  (l-n9)- 
J  {n —  I ):  n! 

We  are  now  in  a  position  to  answer  the  fourth  question. 
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The  probability  that  the  first  empty  interval  of  minimum 
length  6  begins  at  a  distance  I  from  the  origin,  is  equivalent 
to  the  probability  that  the  interval  /  contains  no  void,  while 
the  interval  6  which  immediately  follows  it  is  empty  but 
bounded  on  one  side,  and  this  is 

Fii,  e)j{d)di, 

or,  from  our  definition 

_  8F(l+9,  e) 
61 

From  this  it  follows  that 


dl. 


f. 


F{i,  d)j{e)di  =  F{d,  e)  =  I  -  P{e) 


which  means,  as  was  to  be  expected,  that  the  total  probability 
that  a  void  6  begins  at  a  positive  distance  from  the  origin  is  the 
complement  of  P{0),  which  is  the  probability  that  the  void  6 
begins  precisely  at  the  origin. 

The  probable  delay  in  waiting  for  a  void  of  minimum  length ' 
6  can  therefore  be  written  as 


••CO  ••( 

=  /  iF{i,  d)mdi  =  -  I 


5F(M-M)  ,, 
^ — u — ^^ 


F(/+^,  B)l\    +     /      F(/+^,  Q)dl 

The  first  term  vanishes  at  both  limits,  and  we  also  know  from 
the  properties  of  the  function  F ,  that 

Therefore  we  can  write  the  second  term  in  the  form 

^  F(2g,  Q)  ^  i-P(d) 

^     m       j(d) 

and  this  is  the  answer  to  our  question. 
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It  may  be  noted  that  this  probable  value  p  takes  account  of 
the  cases  where  I  ^  O,  whose  probability  is  P{8). 

We  must  bear  in  mind  that  the  void  6  required  for  a  passing 
maneuver  depends  upon  whether  or  not  the  vehicle  in  question 
has  to  slow  down  before  passing.  Taking  for  example  d^  to  be 
the  void  necessary  for  an  unretarded  passing,  and  60  >  61  to 
be  the  void  required  for  a  retarded  passing,  the  probable  wait- 
ing delay  is  given  by 

^       J  (02)        "'i-Pied 

Fifth  Question:  Mean  number  of  vehicles  contained  in  a  group. 
Probability  of  finding  a  group  with  a  given  number  of  vehicles. 

This  question  is  of  interest  in  the  study  of  the  distribution 
of  spacings.  We  shall  call  a  collection  of  vehicles  a  group  if 
successive  vehicles  contained  in  it  are  separated  by  spacings 
shorter  than  a  given  value  t,  while  successive  groups  are 
separated  by  spacings  longer  than  t.  Now  on  the  basis  of  the 
preceding  developments,  we  know  that  the  probability  of  a 

j(t) 
spacing  longer  than  t  is-  ^  '■,  and  that  the  probable  number  of 

spacings  per  unit  time  is  therefore  J{t)  where  N  is  the  entire 
number  of  spacings. 

This  is  also,  as  we  see,  the  probable  number  of  groups  of 
vehicles  per  unit  time,  and  the  probable  mean  number  of 
vehicles  in  a  group  will  be 

JV 

Moreover,  the  probability  that  there  are  exactly  n  vehicles  in 
a  group  is 

n-l 


a  result  which  remains  valid  for  n  =  1 . 
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These  considerations  afford  a  means  of  determining  the  law 
of  distribution  of  spacings  (by  studying  the  composition  of  the 
various  groups)  of  vehicles  which  pass  a  check  point. 

Conclusion 

It  appears  that  the  results  presented  in  this  paper  can  lead 
to  a  considerably  closer  agreement  between  theory  and  observa- 
tion than  has  been  obtained  heretofore.  Their  interest  consists 
in  the  fact  that  in  spite  of  the  somewhat  complicated  form  of 
the  law  we  have  chosen  for  the  probability  of  spacings  be- 
tween successive  vehicles,  the  values  which  it  gives  in  answer 
to  our  five  fundamental  questions  are  comparatively  simple. 
They  involve  only  the  functions  P,  J ,  and  p.  Of  course  one 
could,  if  one  preferred,  start  with  a  purely  empirical  law,  as 
for  example  the  empirical  curve  of  Figure  2  which  is  simply 

--^,  and  from  this  obtain  the  curves  for  P  and  p  by  graphical 

integration   or  differentiation.    Our  questions  could  thus  be 
answered  without  any  assumed  law  of  probability. 

It  is  hoped  that  these  results  will  be  useful  to  those  interested 
in  traffic  research,  who  recognize  the  growing  importance  of 
statistical  studies  in  this  field.  They  may  also  throw  light  on 
certain  extensions  of  Poisson's  Law. 
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For  readers  interested  in  the  goodness  of  fit  of  the  time-spacing 
distribution  represented  by  the  formula  here  developed,  the 
following  chi-square  (x^)  test  data  are  submitted. 

The  author  has  found  it  desirable  to  replace  the  value  of  the 
various  parameters  given  in  the  text  by  the  following  values: 

/i=i.98  ^2=13.16  7  =  0.583 

e  =  o.8i  1—7  =  0.417 

These  lead  to  the  formula 

(9-0.81)  e 


0.583^  ^-17         +   0.417^  ^3.16 

in  percent  of  the  total  number  of  spacings  (Fig.  1)  and  , 

385  X  io-°-"(»-o-8i)+275Xio-°-°339 

in  number  of  spacings  (table). 

Col.  1  gives  the  limits  of  the  intervals  in  which  the  spacings 
are  classed. 

Col.  2  the  theoretical  number  of  spacings  n,  corresponding  to 
the  entries  in  Column  1 . 

Col.  3  the  theoretical  number  of  spacings  ft  comprised  between 
times  given  by  Column  1 . 

Col.  4  the  observed  number  fo  corresponding  to  ft  of  Column  3. 

Col.  5  the  deviation  ft—fo  • 

Col.  6  the  values  of  -^ —        .  The  sum  of  this  column,  namely 

ft 


I 


Jt 
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FITTING  OF  THE  LAW 

;:   =   385   X    1 0-0 -37  (9-0 -81) -1-2  75  X  I  0-0 -^^Sfl 

(Chi-Square  Test) 


{ft-1o? 

e 

n 

ft 

to 

U-fo 

ft 
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APPENDIX  75 

FITTING  OF  THE  LAW    (Continued) 
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chi-square  =  21.24 
Degrees  of  Freedom  =  25 — 5^20 
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